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We consider the finite clement approximation of the Oldroyd-B system of equations, 
which models a dilute polymeric fluid, in a bounded domain T> C K d , d = 2 or 3, 
subject to no flow boundary conditions. Our schemes are based on approximating the 
pressure and the symmetric conformation tensor by either (a) piecewise constants or 
(b) continuous piecewise linears. In case (a) the velocity field is approximated by con- 
tinuous piecewise quadratics or a reduced version, where the tangential component on 
each simplicial edge (d = 2) or face (d = 3) is linear. In case (b) the velocity field is 
approximated by continuous piecewise quadratics or the mini-element. We show that 
both of these types of schemes satisfy a free energy bound, which involves the logarithm 
of the conformation tensor, without any constraint on the time step f or th e backward 
Euler type time discretization. This extends the results of Boyaval et , ; ; I'IA!0!l on this 
free energy bound. There a piecewise constant approximation of the conformation tensor 
was necessary to treat the advection term in the stress equation, and a restriction on 
the time step, based on the initial data, was required to ensure that the approximation 
to the conformation tensor remained positive definite. Furthermore, for our approxi- 
mation (b) in the presence of an additional dissipative term in the stress equation and 
a cut-off on the conformation tenso r on certain terms in the system, similar to those 
introduced in Barrett and SulJSSSSl for the microscopic-macroscopic FENE model of a 
dilute polymeric fluid, we show (subsequence) convergence, as the spatial and temporal 
discretization parameters tend to zero, towards global-in-time weak solutions of this reg- 
ularized Oldroyd-B system. Hence, we prove existence of global-in-time weak solutions 
to this regularized model. Moreover, in the case d = 2 we carry out this convergence 
in the absence of cut-offs, but with a time step restriction dependent on the spatial dis- 
cretization parameter, and hence show existence of a global-in-time weak solution to the 
Oldroyd-B system with an additional dissipative term in the stress equation. 

Keywords: Oldroyd-B model, entropy, finite element method, convergence analysis, ex- 
istence of weak solutions. 
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2 Jo/m W. Barrett and Sebastien Boyaval 

1. Introduction 

1.1. The standard Oldroyd-B model 

We consider the Oldroyd-B model for a dilute polymeric fluid. The fluid, confined 
to an open bounded domain V C M. d (d = 2 or 3) with a Lipschitz boundary dT>, is 
governed by the following non-dimensionalized system: 

(P) Find u : (t,x) G [0,T) x V i-> as) £ R d , p : (i, as) £ 2? T := (0,T) xDh 
p(t, s)eland(T: (i, as) G [0, T) x 2? (->• cr(t, as) G Rg Xd such that 

Re — + (it • V)tt = - Vp + (1 - e) Am + — div er + f 
\at J Wi 

div u = 

(?(T 1 

— + (u ■ V)<t = (V«)ff + -t(V«) t - — (o- - J) 

«(Q,x) = u°(x) 
cr(0,x) = cr°(x) 
u = 



Here u is the velocity of the fluid, p is the hydrostatic pressure, and er is the 
symmetric conformation tensor of the polymer molecules linked to the symmetric 
polymeric extra-stress tensor r through the relation er = I + ^-t, where 7 is the 
d-dimensional identity tensor and Mg Xd denotes symmetric real d X d matrices. In 
addition, / : (t, x) G 2?t ^ /(£, ^) £ R d is the given density of body forces acting 
on the fluid; and the following given parameters are dimensionless: the Reynolds 
number Re G M>o, the Weissenberg number Wi G K>o, and the elastic-to-viscous 
viscosity fraction e G (0, 1). For the sake of simplicity, we will limit ourselves to the 
no flow boundary condition (ll.lfj) . Finally, Vu(t,x) G R dxd with [Vu]# = ff 1 , 

and (div er)(t, x) G M d with [diver], = £^ =1 ^f. 

Unfortunately, at present there is no proof of existence of global-in-time weak 
solutions to (P) available in the literature. Local-in-time existence results for (P) 
for sufficiently smooth initial data, and global-in-time e xistenc e results for suffi- 
ciently small initial data can be found in Guillopc and 

SauiGSSQ] 

for a Hilbert space 

framework, and in Fernandez- Cara et a/ F ^ for a more general Banach space 

framework. Global-in-time existence results for the corotational version of (P); that 
is, where Vu in (|l.lc[) is replaced by its anti-symmetric part i(Vu — (Vit) T ) can 
be found in Lions and MasmoudiP^^ We note that such a simple change to the 
model leads to a vast simplification mathematically, but, of course, it is not justi- 
fied on physical grounds. Finally, global- in-time existence results for (P) in the c ase 
/ = and for initial data close to equilibrium can be found in Lei et a /j LLZ08 ] 

This paper considers some finite element approximations of the Oldroyd-B sys- 
tem, possibly with some regularization. In the regularized case, we show (subse- 
quence) convergence of the approximation, as the spatial and temporal discretiza- 
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Existence and Approximation of a (Regularized) Oldroyd-B model 3 

tion parameters tend to zero, and so establish the existence of global-in-time weak 
solutions of these regularized versions of the Oldroyd-B system. The first of these 
regularized problems is (P Q ) obtained by adding the dissipative term a Act for a 
given a £ M>o to the r ight-ha nd side of (|l.lc|) . as considered computationally in 
Sureshkumar and BerispSSH w ith an additional no flux boundary condition for 
cr on dV. The second is (P„) where, in addition to the regularization in (P a ), 
the conformation tensor cr is replaced by the cut-off f3 L (cr) on the right-hand side 
of (|l.la|) and in the terms involving u in (jl.lcj) . where f3 L (s) := min{s, L] for a 
given L » 1. Similar regularizations have been introduced for the microscopic- 
macroscopic dumbbell model of dilute polymers with a finitely extensible nonlinear 
elastic (FENE) spring law, see Barrett and SulipHflS an d f or the c onverg ence of 
the finite element approximation of su ch models, see Barr ett and Suli.ESlQI i n f act 
it is argued in Barrett and Siil -BS07 ] an( ^ SchiebeJ ^ c ^^ l that the dissipative term 



a Act is not a regularization, but is present in the original model with a positive 
a -C 1- Here we recall that the Oldroyd-B system is the macroscopic closure of 
the microscopic- macro scopic dumbbell model with a Hookean spring law, see e.g. 
Barrett and SuliPSEl 

Overall the aims of t his paper are threefold. First, we extend previous results 
in Boyaval et a ;| BLM0 9] f or a g n jt e element approximation of (P) using essentially 
the backward Eulcr scheme in time and based on approximating the pressure and 
the symmetric conformation tensor by piecewise constants; and the velocity field 
with continuous piecewise quadratics or a reduced version, where the tangential 
component on each simplicial edge (d — 2) or face (d — 3) is linear. We show 
that solutions of this numerical scheme satisfy a discrete free energy bound, which 
involves the logarithm of the conformation tensor, without any constraint on the time 
step, whereas a time constraint based on the initial data was required in Boyaval 
et a /j BLM0 9l m orc [ er i ensure that the approxima tion to the conformation tensor 
cr remained positive definite. See also Lee and XupHU w here the difficulties of 
maintaining the positive definiteness of approximations to cr are also discussed. 
We achieve our result by first introducing problem (Pg), based on a regularization 
parameter S <E K.>o. (Pa) satisfies a regularized free energy estimate based on a 
regularization of In and is valid without the positive definiteness constraint on the 
deformation tensor. 

Second, we show that it is possible to approximate (P) with a continuous (piece- 
wise linear) approximation of the conformation tensor, such that a discrete free 
energy bound still holds. We note that a piecewise constant approximation of the 
conformation tensor was necessary in Boyaval et Q ;j BLM0 9] j n orc j er ^ treat the 
advection term in (|l.lc|) and still obtain a discrete free energy bound. 

Third, we show (subsequence) convergence, as the spatial and temporal dis- 
cretization parameters tend to zero, of this latter approximation in the presence 
of the regularization terms stated above to global-in-time weak solutions of the 
corresponding regularized form of (P). 

The outline of this paper is as follows. We end Section [T] by introducing our 
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4 John W. Barrett and Sebastien Boyaval 

notation and auxiliary results. In Section [2] we introduce our regularizations of In 
based on the parameter S e (0, |] and the cut-off L > 2. We introduce our regu- 
larized problem (Ps), and show a formal free energy estimate for it. In Section [3j 
on assuming that I? is a polytope, we introduce our finite element approximation 
of (Ps), (P^fc) based on approximating the pressure and the symmetric conforma- 
tion tensor by piecewise constants; and the velocity field with continuous piecewise 
quadratics or a reduced version, where the tangential component on each simplicial 
edge (d = 2) or face (d = 3) is linear. Using the Brouwer fixed point theorem, we 
prove existence of a solution to (Pf'fy and show that it satisfies a discrete regular- 
ized free energy estimate for any choice of time step; see Theorem 13. II We conclude 
by showing that, in the limit 6 — > + , these solutions of (Pfl) converge to a so- 
lution of (P^ ) with the approximation of the conformation tensor being positive 
definite. Moreover, this solution of (P^*) satisfies a discrete free energy estimate; 
see Theorem [32] 

In Section 2] we introduce our regularizations (P^) of (P) involving the dissi- 
pative term a Aer on the right-hand side of (|l.lc[) . and possibly the cut-off /3 L (cr) 
on certain terms involving a in (JTTTaJc). We then introduce the corresponding 
regularized version (Pj^j), and show a formal free energy estimate for it. In Sec- 
tion [5] we introduce our finite element approximation of (P^j), (^aSh*) based on 
approximating the pressure and the symmetric conformation tensor by continuous 
piecewise linears; and the velocity field with continuous piecewise quadratics or the 
mini-element. Here we assume that the finite element mesh consists of non-obtuse 
simplices. Using the Brouwer fixed point theorem, we prove existence of a solution 
to (P„ a ft *) and show that is satisfies a discrete regularized free energy estimate for 
any choice of time step; see Theorem l5.ll We conclude by showing that, in the limit 
6 — > + , these solutions of (P^/f*) converge to a solution of (P^ A *) with the 
approximation of the conformation tensor being positive definite. Moreover, this 
solution of (P^A ) satisfies a discrete free energy estimate; see Theorem 15.21 

In Section[6]we assume, in addition, that I? is a convex polytope and that the fi- 
nite element mesh consists of quasi-uniform simplices. We then prove (subsequence) 
convergence of the solutions of (P^^ 4 ), as the spatial and temporal discretization 
parameters tend to zero, to global-in-time weak solutions of (P„); see Theorem l6.2l 
Finally in Section [3 on further assuming that d = 2 and a time step restriction 
dependent on the spatial discretization parameter, we prove (subsequence) conver- 
gence of the solutions of (P^), as the spatial and temporal discretization param- 
eters tend to zero, to global- in-time weak solutions of (P Q ); see Theorem 17.21 We 
note that these existence results for (PqP^) are new to the literature. In addition, 
the corresponding L°°(0, T; L 2 (fl)) n L 2 (0, T; H 1 ^)) norms of the velocity solution 
of (PqP^) are indepen dent of the regularization parameters a (and L). 

In a forthcoming paper JEBQSI we w in extend the ideas in thi s pape r to a related 
macroscopic model, the FENE-P model; see Hu and LelievreJSG^I w here a free 
energy estimate is developed for such a model, as well as Oldroyd-B. In addition, 
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Existence and Approximation of a (Regularized) Oldroyd-B model 5 

we will report in the near future on numerical computations based on the finite 
element approximations in this paper and Barrett and Boyaval|EEQ9| 

1.2. Notation and auxiliary results 

The absolute value and the negative part of a real number s £ K are denoted by 

If 

Vi 4-Vid o/a-f cirin tyi cif fin nno 

^SPD 



\s\ := max{s, — s} and [s]_ = min{s, 0}, respectively. In addition to R^ xd , the set 
of symmetric W dxd matrices, we let 1™ be the set of symmetric positive definite 



T>dxd 



matrices. We adopt the following notation for inner products: 

d 

v ■ w :— ViW i = v T w — w T v Vi>, w £ R d , (1.2a) 

i=l 
d d 

<t> * : = E E tutu = tr (^) = tr (^) v< ^ * e Rdxd > ( L2b ) 

«=1 3=1 

Vo:: V<..-: EE Vo -' V ' ; ''/ V0,0£R dxd ; (1.2c) 

«=1 3=1 

where - T and tr (■) denote transposition and trace, respectively. The corresponding 
norms are 

|M| :=(«•«)*, ||Vt>|| := (Vw : V©)* Vw £ R d ; (1.3a) 
||0|| := (0:0)4, ||V0|| :=(V0:: V0)s, V0 £ K dxd . (1.3b) 

We will use on several occasions that tr(0) = tr(0 T ) and tr(00) = tr(00) for all 
0,0 £ R dxd . In particular, we note that: 

0X T :0 = X0:0 = x:00 V0, £ M| x d , X £ M dxd , (1.4a) 
H00II < ||0|| ||0|| V0,0£R dxd . (1.4b) 

In addition, for any € M^ x<i , there exists a diagonal decomposition 

= O T DO => tr(0)=tr(D), (1.5) 

where O £ R dxd is an orthogonal matrix and D £ K dxd a diagonal matrix. Hence, 
for any g : R ->■ M, one can define g(0) £ Rg Xd as 

g(0) := T 3 (D)0 =► tr ( 3 (0)) = tr ( 3 (D)) , (1.6) 

where 5(D) £ Rg Xd is the diagonal matrix with entries [<?(D)]jj = g(Du), i = 1 d. 
Although the diagonal decomposition (jl.5[) is not unique, (jl.6l) uniquely defines 
g{<$>). Similarly, one can define <?(0) £ R^, when £ and <? : M >0 -> R. 

We note for later purposes that the choice g(s) = \s\ for s £ R yields that 

d _1 (tr(|0|)) 2 < ||0|| 2 < (tr(|0|)) 2 V0 £ Mf d . (1.7) 
We adopt the standard notation for Sobolev spaces, e.g. i? 1 (2?) := {77 : V — > 
R : /dII^I 2 + ll V? ?ll 2 ] < °°} with H o( v ) bein S the closure of C£°(£>) for the 
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6 John W. Barrett and Sebastien Boyaval 

corresponding norm || • ||jyi(x>)- We denote the associated semi- norm as | • \h 1 {t>)- 
The topological dual of the Hilbert space Hq(V), with pivot space L 2 (T>), will be 
denoted by 7J _1 (2?). We denote the duality pairing between i? -1 (X>) and Hq(T>) 
as (•, -)_ffi(p)- Such function spaces are naturally extended when the range M is 
replaced by R d , R dxd and WL dxd ; e.g. H 1 ^) becomes [H 1 (V)] d , [H 1 (V)] dxd and 
[H 1 {V)] d s xd , respectively. For ease of notation, we write the corresponding norms 
and semi-norms as || • ||iji(x>) and \-\h 1 (v)i respectively, as opposed to e.g. || ■ || [H 1 (T>)] d 
and | • \[H 1 (v)] d j respectively. Similarly, we write (-,-)h^(V) f° r t ne duality pairing 
between e.g. [H~ 1 {V)] d and [HliV)] d . We recall the Poincare inequality 

/ \\v\\ 2 <C P I ||Vu|| 2 V«e[ff W, (1-8) 

JV JT> 

where Cp G R>o depends only on T>. For notational convenience, we introduce 
also convex sets such as [H 1 (V)]f p d D := {(f) e [H 1 {V)] dxd : e R d s x p d D a.e. in V}. 
Moreover, in order to analyse (P), we adopt the notation 

W := [Hl{V)] d , Q:=L 2 (V), V := jw E W : J q div v = V</ G q| , 

S := [LHV)]** d and S PD := [L\V)] dx p d D . (1.9) 

Finally, throughout the paper C will denote a generic positive constant independent 
of the regularization parameters 5, L and a; and the mesh parameters h and Ai. 



2. Formal free energy estimates for a regularized problem (Ps) 
2.1. Some regularizations 

Let G : s € M>o H Ins S R be the logarithm function, whose domain of definition 
can be straightforwardly extended to the set of symmetric positive definite matrices 
using (|1.5|) . We define the following two concave C 1 (R) regularizations of G based 
on given parameters L > 1 > 5 > 0: 

(G(s), Vs><5 , „ T \g l (s), Vs>6 

\§ + G((5)-l, Vs<5 \G s (s), Vs<S 



T \f + G(L) - 1 , Vs > L 

where G L : s £ R >0 ^ < L 

\G(s), Vse(0,L] 

We define also the following scalar functions 
(3l L \ s ) := (g< L) '(s)) _1 Vs e R and /3 (i) (s) := (g (l) '(s)) 



(2.1) 



Vs e K >0 ; 

(2.2) 



where, here and throughout this paper, -W denotes an expression with or without 
the superscript and a similar convention with subscripts. Hence we have that 

ft:seM4 max{s, 6} , ft 1 : s 6 1 4 min{^(s), L} , 
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Figure 1. The function G and its regularizations. 



P : 8 G 



and 



/3 L : s G M>o H- min{/3(s), L} . 



We note for example that 

|| || 2 < dL 2 V<p e R d / d and ||^ L (</>)|| 2 < dL 2 V0 G Rf P d D • 
Introducing the concave G 1 (R) functions 

H^(s) :=G 5 L -\{s) VsGM and H s (s) := G 8 ' 1 (s) Vs G 
it follows from (j2~Tj) and that 

fff ) '(Gf ) '( S ))=4 i) ( S ) feel. 
For later purposes, we prove the following results concerning these functions. 

Lemma 2.1. The following hold for any <t>, i/> G R d s xd and for any L > 1 > S > 
that 



(2.3) 
(2.4) 
(2.5) 
(2.6) 



tr 



[^ ) (0)][G^ ) '(0)] = [Gf ) '(0)][ y 8f ) (0)] 
{^\<j>) + \!3f\ct>)]- l -2l) >0, 



trU-Gf^)-/) >0, 



0-/3< L) (0)) : (l-G< L) '(0)) >0, 

(0 - V) : (flf >'(tf)) > tr (g< L) (</>) - G< L) (</>)) , 
(0 - V) : (g1 L) '(0) - Gi L) '(^)) > 5 2 |G< L) '(0) - G< L) '(</>) 



(2.7a) 
(2.7b) 
(2.7c) 
(2.7d) 
(2.7e) 
(2.7f) 
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8 John W. Barrett and Sebastien Boyaval 

In addition, if 5 £ (0, |] and L > 2 we have that 
tr 



U-Gf\<t>)) >{^\ and 0:(j-G< L) '(0))>§|M|-d. (2-8) 



Proof. The result (|2.7a|) follows immediately from (jl.6[) and as j3g L ^ (s) = G S L ^ (s) 
for all s€l. The desired results |277b]-d) follow similarly, on noting the scalar in- 
equalities pf ] (s) + [pf ] (s)]" 1 > 2, s-Gf\s) > 1 and {a-pf 3 (s))(l-Gf 3 '(s)) > 
for all s G R. 

We note that G S L ^ are concave functions like G, and hence they satisfy the 
following inequality 

(si-s 2 )G g L) '(s 2 ) >G s L \ Sl )-G g L) {s 2 ) V Sl ,s 2 GM; (2.9) 

Hence for any 0, ip G Rg Xd with = O^D^O^ and -0 = O^D^O^,, where 
0^,0,/, G R dx<J orthogonal and D ,D^ G R dxd diagonal, we have, on noting the 
properties of trace, that 

(0 - 0) : G^'(0) = tr ((0 - 0)Gf >'(0)) = tr ((O T D O - D^Gf } '(D^) 

(2.10) 

where O = O O^ G R dxd is orthogonal and hence E^=i[°y] 2 = Etil ^] 2 = 1 
for j = 1 —> d. Therefore we have, on noting these properties of O, (|2.9p and (|1.6p . 
that 

did \ 

tr ((O T D O - D^)Gf > (D*)) = £ ^[O,,] 2 ^],, - [D^ [G^'(D^)]« 

»=i y=i / 

d d 

»=i j=i 

>EEN 2 (^ L) (D )] w - [G< L) (D^) 

i=l 3=1 



tr (c^D^-tr (Gf(D^)) 

(Gf>(0)-G< L >(0)) • (2- 11 ) 



= tr 

Combining ((2~TU|) and (ETTTI) yields the desired result (l2T7cl) . 

We note that —G S L ^ G G 0,1 (R) is monotonically increasing with Lipschitz con- 
stant (5~ 2 and so 

- s 2 )(G s Ly ( Sl ) - G s L) \s 2 )) > 6 2 [G 5 L) '( Sl ) - G 5 L) '(s 2 )} 2 V Sl , s 2 e R. 

(2-12) 

Then, similarly to (|2~10l) and (|2~TT|) . we have, on noting (f2TT2]) . that 
-(0-0):(G< L) '(0)-G< L) '(0)) 
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[tr ((D - OD^O T )G[ L) '(D^)) - tr ((O T D O - D^G^D 



= -EEW 2 ~ (lG<s L) '(n<t>)h ~ [G^'(D^)]«) 

»=i j=i 

i=i j=i 

= S 2 tr ((G< L) '(0) - Gf y W)f) = 5 2 \\G ( s Ly (4>) - G< L) M|| 2 (2-13) 

and hence the desired result (|2.7f[) . 

Finally the results (12.81) follow from (| 1 . 6[) and (|1.7p on noting the following scalar 
inequalities 

s-Gf 3 (s) > ^ 2 and s(l-G^ 5 (s)) > Vs6l, (2.14) 

which are easily deduced if S € (0, |] and L > 2. □ 

Clearly (|2.7e|) holds for any concave function g g G 1 (M), not just G^ , and this 
implies that 

(0-V'):ff / W>tr(ff(0)-5(^))>(^-^):fl / (^) V0,</>eMf d . (2.15) 

For a convex function 5 e G 1 (R), the inequalities in (|2.15[) are reversed. Hence for 
any concave or convex function g e G 1 (M) and for any 4> S G 1 ([0, T); M^ xrf ) one 
can deduce from the above that 

d - tr (g(d>)) = tr ( ^-g'(<P)) = </(<£) Vt £ [0, T). (2.16) 



\ d< J \ dt J 

Of course, a similar result holds true for spatial derivatives. Furthermore, these 
results hold true if <p is in addition positive definite, and g € G 1 (K>o) is a concave 
or convex function. Finally, we note that one can use the approach in (|2.1ip to 
show that if g £ G 0,1 (R) with Lipschitz constant gup, then 

\\g(4>)-gm\<gu P \\4>-ip\\ v<f>^eRf d . (2.17) 

2.2. Regularized problem (Ps) 

Using the regularizations G$ introduced above with parameter S we consider the 
following regularization of (P) for a given 5 £ (0, i]: 

(Ps) Find us : (t,x) € [0,T) x V ^ u s (t,x) e R d , ps : (t,x) e(0,T)xI)4 
p 5 {t, x) £ R and cr s : (t, x) e [0, T) x 2? h> <r 5 (t, x) 6 M^ xd such that 

Re + (u s ■ V)u s J = -Vps + (1 - e)Au s + divft(<r tf ) + f 

on V T , (2.18a) 
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10 John W. Barrett and Sebastien Boyaval 

div« 5 = onX> T , (2.18b) 

+ (u 5 ■ V)a s = (Vti^^j) + Ps{ct s ){Vus) t - — (a s - I) 

on V T , (2.18c) 

u 4 (0, x) = u°(x) VxeV, (2.18d) 

0-5(0, x) = a-°(x) \fx G T> . (2.18e) 

u s = on (0,T) x dV. (2.18f) 

2.3. Formal energy estimates for (Ps) 

In this section, we derive formal energy estimates, see e.g. (|2.21[) below, where we 
will assume that the triple (us,ps,as), which is a solution to problem (Ps), has 
sufficient regularity for all the subsequent manipulations. 
We will assume throughout that 

/ e L 2 (0,T; [H-\V)] d ) , u° £ [L 2 (V)] d , and <x° E [L°°(V)] d % 

with a° min \\$\\ 2 <£ T <r°(x)£ << ax IIIH 2 V|eM rf for a.e. x in V; (2.19) 

where < in , < ax £ R >0 . 

Let i^(it5,<r 4 ) denote the free energy of the solution (us,ps,o~s) to problem 
(P 4 ), where F s (; ■) : W x S -)• R is defined as 

F,(«, 0) := ^ jf H| 2 + ^ jf tr(0 - G 4 (0) - I) . (2.20) 

Here the first term ^ ||f|| 2 corresponds to the usual kinetic energy term, and 
the second term J v tr(<p — Gs(4>) — I) is a regularized version of the relative 
entropy term tr( </> — G(<p) — I) introduced in Hu and LelievrepEEl see a i so 

Jourdain et ^ fjLLO06 | 

Proposition 2.1. Let (us,Ps, c 4 ) be a sufficiently smooth solution to problem (Ps)- 
Then the free energy Fs(us, as) satisfies for a.e. t G (0, T) 

yFs(us,<rs) + (l-e) [ \\Vu 5 \\ 2 + -^[ tr(/3 4 (o- 4 ) + tfsfa)]' 1 - 21) 
at J v 2Wr Jv 



v 

< (f, us) Hl(p)- (2.21) 



Proof. Multiplying the Navier-Stokes equation with us and the stress equation 
with 2^(1 — G' s (<Ts)), summing and integrating over T> yields, after using integra- 
tions by parts and the incompressibility property in the standard way, that 

5J!-||m 4 || 2 + (1 - £)||Vm 4 || 2 + 4 ? Ps(as) : Vu 5 



2 dt" " v /n " Wi 



2Wi 



d \ 1 

—as + (us ■ V)a s J + — (a 5 - 1) 



(I-G' s (as)) 



1 



1 
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- J [(Vu s )(3s(cts) + Ps^s)(Vus) T ) : (I - G' s (a s )) = (f,u s ) Hl{v) . 

(2.22) 

Using (|2.16l) and its spatial counterpart, we first note that 

^as + (u s ■ V)<r a ^ : (I - G' s (a s )) = (J^ + (u s ■ V)j \x(tr s -G s (*s)) ■ 

(2.23) 

On integrating over T>, the (us ■ V) part of this term vanishes as u§(t, •) G V. On 
using trace properties, (|2.7a[) and the incompressibility property, we obtain that 

((Vu s ) p B (trs)) ■ (I - G' 5 (<ts)) = tr ((Vu s ) (3 5 (a 5 ) - (Vu s ) Ps(cts)G's(<ts)) , 

= tr((V« s )/3 s (<7 4 ) - Vu s ) , 
= tr ((Vus) fis(crs)) ~ div us , 
= tr ((Vu s ) j3 5 (a s )) ■ (2.24) 
On noting (|1.4a[) . we have also that 

(Psfa) (Vusf) : (I - G' s (as)) = tr ((Vu s ) Ps(rr S )) ■ (2.25) 

Therefore the terms involving the left-hand sides of (|2.24l) and (|2.25p in (12.221) cancel 
with the term ^f3s(crs) ■ Vus in (|2.22l) arising from the Navier-Stokes equation. 
Finally, for the remaining term we have on noting (ll.2b|) and (|2.7a[ d) that 

(a s -I): (I- G's(crs)) = \{a s - Ps{crs)) + (Ps(^s) - I)] ■ (I - G' s (a s )) 
>(Ps(<ts)-I):(I-G' 5 (<ts)) 

= tr(p 5 ((Ts) + McTsT 1 - 21) . (2.26) 
Hence we obtain the desired free energy inequality (|2.21[) . □ 

Corollary 2.1. Let (us,ps, &s) be a sufficiently smooth solution to problem (P5). 
Then it follows that 

sup F s (us(t, -),(Ts(t, ■)) 

te(o,T) 

\(l - e)\\Vu 5 f + ^tr(Ps(crs) + M'tT 1 ~ 2I ) 



< 2 (f s (u°, <x°) + 1±^L \\ff WiB - lm 



(2.27) 



Proof. Smooth solutions (us,ps, crs) of (P5) satisfy the free energy estimate (|2.20j) . 
One can bound the term (/, us) hI(t>) there, using the Cauchy-Schwarz and Young 
inequalities for v £ M>o, and the Poincare inequality (|1.8[) . by 



1 2 2 
(f,us) H i(v) < ll/llff-i(x>) \\ u s\\m(v) ^ ^2 Wf\\H-i(v) + Y W Us \\hhv) 



1 



1 
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12 John W. Barrett and Sebastien Boyaval 

< ^ \\ff H - Hv) + y(l + C P ) \\Vu & \\l 2{v) . (2.28) 

Combining (I2.28|) and (|2.20l) with v 2 = (1 — s)/(l + C P ), and integrating in time 
yields the desired result (|2.27[) . □ 

We note that the right-hand side of (|2.27p is independent of the regularization 
parameter 6 if er° is positive definite. 



3. Finite element approximation of (Ps) and (P) 
3.1. Finite element discretization 

We now introduce a finite element discretization of the problem (Ps), which satisfies 
a discrete analogue of (|2.2ip . 

The time interval [0,T) is split into intervals [f 1 - 1 ,t n ) with At n = t n - f -1 , 

n = 1, . . . ,Nt- We set At := max„ = i n t At n . We will assume throughout that 

the domain T> is a polytope. We define a regular family of meshes {Th}h>o with 
discretization parameter h > 0, which is built from partitionings of the domain T> 
into regular open simplices so that 

— N K hk 

T> = Th '•= U Kk with max — < C . 

k=l k=l....,N K Pk 

Here pk is the diameter of the largest inscribed ball contained in the simplex Kk 

and hk is the diameter of Kk, so that h = maxfc = i n k hk- For each element 

Kk, k = 1, . . . ,Nk, of the mesh Th let {-Pj fc }f = o denotes its vertices, and {n^}f =0 
the outward unit normals of the edges (d = 2) or faces (d — 3) with n\ being 
that of the edge/face opposite vertex P/% % = 0, . . . , d- In addition, let {Vi{ x )}i=Q 
denote the barycentric coordinates of x 6 Kk with respect to the vertices {P^ }f—o', 
that is, r\\ E Pi and 77* (P*) = 5y, i, j = 0, . . . , d. Here P m denote polynomials of 
maximal degree m in x, and <5y the Kronecker delta notation. Finally, we introduce 
dTh '■= {Ej}f=i as the set of internal edges Ej of triangles in the mesh Th when 
d = 2, or the set of internal faces Ej of tetrahedra when d = 3. 

We approximate the problem (Ps) by the problem (P^h) based on the finite 
element spaces W° x Q° x S°. As is standard, we require the discrete velocity- 
pressure spaces W° x Q" C W x Q satisfy the discrete Ladyshenskaya-Babuska- 
Brczzi (LBB) inf-sup condition 



q divv 

inf sup — — |j — ji > p± > , (3.1) 

9GQ°uew° Il9lli2(p) II u IIhi(x>) 

see e.g. pll4 in Girault and Raviart JSEHS] j n following^ we set 

:= W^ C W or W^'~ C W , (3.2a) 

Q° :={?eQ : ?|jf h GP fc=l,...,JVjc}cQ, (3.2b) 



1 



1 
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and S° h :={4>eS : cf> \ Kk E [Po]f d k = 1, . . . , N K } C S; (3.2c) 

where 

wJ-{«e[C(P)fnw :() | Kl e[P 2 f k = l,...,N K }, (3.3a) 

W£ 1- := G [C(X>)] d nWiw^G [P^ © span{<rf }f =0 fc = 1, . . . , JV K } • 

(3.3b) 

Here, for k = 1, . . . , A^- and i = 0, . . . , d 

d 

<i k i {x)=n k i J] t$(x) foYXEK k . (3.4) 

We introduce also 



which approximates V. I t is we ll-known the choices (13.2a[ b) satisfies (|3. 1|1 . see e.g. 
p221 in Brezzi and Fortir l BF92 l for W° h = W%, a nd Chapter II, Sections 2.1 (d = 2) 
and 2.3 (d = 3) in Girault and RaviartPEESI for w o = w£' _ . Moreover, these 
particular choices of S° and Q° have the desirable property that 
<PeS° h => I-G' s (d>) eS° h and tr(0-G 5 (0) -I) G Q° h , (3.5) 

which makes it a straightforward matter to mimic the free energy inequality (|2.2ip 
at a discrete level. Since S^ is discontinuous, we will use the discontinuous Galerkin 
method to approximate the advection term (u$ ■ V)c,5 in the following. Then, for 
the boundary integral s, we w ill make use of the following definitions (see e.g. p267 
in Ern and Guermone ffi^^ ) . Given v G W°, then for any cp G S° (or Q°) and for 
any point x that is in the interior of some Ej G dl~h, we define the downstream and 
upstream values of d> at x by 

4> +v {x) — lim <p(x + pv(x)) and <fi~ v (x)= lim <p(x + pv(x)) ; (3.6) 

p->0+ p->0~ 

respectively. In addition, we denote by 

ld>Uv(x)=d> +V (x) -<!>-* \x) and {<t>Y{x)= <i>+V{x)+ 2 Ct) (3.7) 

the jump and mean value, respectively, of <f> at the point x of boundary Ej. From 
(|3.6[) . it is clear that the values of 4> +v \E j and 4>~ v \E j can change along Ej G <977i. 
Finally, it is easily deduced that 

We . N K 



j2 / i« • ni[«ii-M, <?r = - £ / * ^ 

j=i fc=i 



W G W£, q u q 2 G Q£ ; (3.8) 

where n = n(Ej) is a unit normal to Ej, whose sign is of no importance, and TiK k 
is the outward unit normal vector of boundary dK^ of Kk- We note that similar 
ideas a ppear in upwind schemes; e.g. see Chapter IV, Section 5 in Girault and 



appear 
rt GR86 



Raviao ^^"" ! for the Navier-Stokes equations. 



i 



i 
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3.2. A free energy preserving approximation, (Pg^), of (Ps) 

For any source term f E L 2 (0,T; [H- 1 (V)] d ) , we define the following piecewise 
constant function with respect to the time variable 

rt" 



/ At ' + (V) = /"(•) := J f(t,-)dt, t 6 [t n -\t n ), 



It is easily deduced that 



n = l,...,N T . 



(3.9) 



and 



X>i n ||rilff-i (2 >) < / \\f(t,-W H - H v)dt for any re [1,2] 
At ' + -> / strongly in L 2 (0,T; [i2'- 1 (Z>)] d ) as At 0+ . 



(3.10a) 
(3.10b) 



Throughout this section we choose tt° £ V° to be a suitable approximation of 
it such as the L 2 projection of u° onto . We will also choose cr° g to be the 
L 2 projection of <r° onto S° . Hence for As = 1, ... , -/Vr- 



1 



(3.11a) 



\ K k\ J Kk 

where \Ku\ is the measure of Kk; and it immediately follows from (|2.19p that 

<in U\\ 2 < £ T °l \K h £ < UW 2 ^eR d . (3.11b) 

Our approximation (Pfl) of (Pa) is then: 

(Pfl) Setting = K,<r°) G V° x S°, then for n = 1,...,7V T find 

(«j fc , o-^) G V° x S° such that for any test functions (v,tj>) G V° x S° 



.1—1 " 



Re 



Re 



(K^wJ-t,-ttj fc .((u^.v)»)' 

(r,«) ff i(i,), (3.12a) 



(1 - e)V«?, h : V« + Wi S (<rlh) : 



At n 



E 



0-2((V<J/3 5 (< /l )) :</)- 



- 1) ■ <i> 







Ie, 









(3.12b) 



Wi 
=0 . 

In deriving (Pfl), we have noted (|1.4aj) and that 

/ v[(z-V)w}=- I w-[(z-V)v] VzgV, Vv,w G [H l (D)] d . (3.13) 
Jv Jv 



our 



Once again we refer to p267 in Ern and 

GuermoncJESE] 

for the consisten cy of 
stated approximation of the stress convection term, see also Boyaval et a /| BLM0 9] 



i 



i 
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Before proving existence of a solution to (P^) , we first derive a discrete analogue 
of the energy estimate (|2.21[) for (P$£); which uses the elementary equality 



2si(si - s 2 ) = s 2 — s 2 + (si — s 2 ) 2 Vsi, s 2 G M. 



(3.14) 



3.3. Energy bound for (Pf^) 

Proposition 3.1. For n = 1,...,Nt, a solution (ug th , tr^jA G V° x 
\3.12a\ b ), if it exists, satisfies 



to 



A/„ ' 2Ai„V' " AVl 



It"- 1112 



(1-e) / ||V< ; , 
v 



n 1 1 2 



+ 



2Wi 2 



/ tr(/3,( < r^) + [^«,)]- 1 -2/) 



< (f\ui h ) Hh(v) < \o-e)J iiv^f+i^nrii^ 



i(x>) ■ 



(3.15) 



Proof. Similarly to the proof of Proposition 12.11 we choose as test functions v = 
ul h G V° h and 4> = ^ (i - Ggi^lfS) e S ° in <fS32al b). and obtain, on noting 
dOH) and ([2~7al d). that 

(f"i u th)Hi(V) 



> 



V 



Re /" KJP^K^f 
2 I Ai„ 



+ 



2Wi 

e 



At r 



(J-G^)) 



(1- £ )|]V< 112 



b S,h\ 



. tr(/3 5 « /l ) + [/3 5 « /l )]- 1 -2/) 

Z Wl Jx> 



N E 

—y 



(J - Gi« h )) 



(3.16) 



We consequently obtain from (|3.16j) . on noting (jl.2b|) and (|2.7c[) applied to the edge 
terms as well as the discrete time derivative term for the stress variable, that 



</":'"5,/ l )ff 1 (23) 

> 



Re /K^-KiT 



At r 



Kh ^-)+il-e)\\V»U\' 



Atr, 



+ 



+ 



tr 



2Wi Jv 

e 

2Wi 2 j r 



(*Ih - Gs(< h )) tr (og - g^rg)) 
Ai„ 



tr(/3 5 « h ) + [^« h )]- 1 -2J) 



i 



i 
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16 John W. Barrett and Sebastien Boyaval 



—t 



3=1 " E i 



[tr «,-G 5 «J 



n-1 



(3.17) 



Finally, we note from (l3~8")l . ([3"T5]l and asu^ eVj that 



N E 

E 



[tr (al h -G s {al h 

N K 



k=l hK k v ' 

AT* „ 

-£/ divfti^tr^-G,^)) 
- / tr^-Gjf^divu^^O. 



(3.18) 



Combining p,17[) and (|3.18p yields the first desired inequality in Q3.15P . The second 
inequality in p,15[) follows immediately from (|2.28[) with v 2 = (1 - e)/(l + Cp). □ 

3.4. Existence of a solution to (P^) 

Proposition 3.2. Given (w^ 1 ,^ 1 ) eVjxSj and for any time step At n > 0, 
then there exists at least one solution {u^h^lhj ^ x t° I3.12al b). 

Proof. We introduce the following inner product on the Hilbert space V° x 

((to, </>), (v, cj>)) v = f [wv + il>:ct>] V(t», V), (», 0) e V?, x S£ . (3.19) 

Given (uj^ 1 ,^ 1 ) 6 V° x S°, let J : V» x S» 4 V° x S° be such that for any 
(w^)eV h xS° h 



Re 



it) — it 



i-l ' 



Re r 



- (1 - e)Vw : Vv + JLp s (i/,) : Vv 



2((V™)/? 5 (t/O):0 + — (</>-/):</> 
Wi 

N E 

E 



:0 



V(u, 0) € V" x S° . (3.20) 



i 



i 
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We note that a solution (ug h , (Tg h ) to (|3.12a[ bV il it exists, corresponds to a zero 
of J 7 ; that is, 



{F{ul h ,al h ),{v,4>)) 



v 







V(«,0) G V° h x Si 



(3.21) 



In addition, it is easily deduced that the mapping T is continuous. 

For any (w,ip) G V° x S° , on choosing (i>, 0) = (to, (I ~ G' s (ip)) , we obtain 
analogously to (|3.15|) that 



(^(™,t/0,(™,^(J-G' 5 (V0) 



> 



2Wi 

F 4 (to,^)-F a (u 



n— 1 n— 1 



Re 



.n-l il 2 



1-E 



2Wi 2 



At n 2At n J v 



\Vw\ 



'irWl-^v) ■ (3-22) 



2(1 -e) 

Let us now assume that for any 7 G K>o, the continuous mapping J- has no zero 
(ug h , (TgfJ satisfying (|3 . 2 1 1) . which lies in the ball 



B 1 :={(v,<j>)eV° h xSl : \\(v, 0)|| D < 7} I 



where 



\(v,(f>)\\ D := [((v,<l>),(v,<l>)) D ]* = 



\\4>\\ 



(3.23) 



(3.24) 



Then for such 7, we can define the continuous mapping Q~ : B 1 — > B 1 such that for 
all (v,4>) G B 1 



mv,<t>)\\v 



(3.25) 



By the Brouwer fixed point theorem, C? 7 has at least one fixed point (u? 7 ,V' 7 ) in 
B T Hence it satisfies 

||(>7><Mlx> = \\^i W l^l)\\T> = T 

It follows, on noting l[5^c| . (|3T2"4l and ([335]) . that 
1 



(3.26) 



7lli°°(X>) 



< 



min^jvj, \K k \ J v 



IIV-, 



(3.27) 



where /i h := [l/(min feeA r K |-f£T fe |)]i It follows from ([2~^0jl . (|2T5]l . (f3T27|) and ([3226} 
that 



2Wi 



> 



> 



Re 

T 

Re 



+ 



4Wi 



4Wi/i ft 7 



tr(0 7 - G a (0 7 ) - I) 
||^ 7 ||-2d|2?| 
ll^ 7 IU~(i>) / I |V'- 



v 



ed\V\ 
2Wi 



1 



1 
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,'Re e 
> mm — , 



/[ 



2 mifihjj \Jv 



Re 



ed\V\ 



2Wi 



sd\V\ 
2Wi 



2 '4Wi/i /l7/ 

Hence for all 7 sufficiently large, it follows from (|3.22p and p.28[) that 
(T(w^), (w„ ^ (J- G> 5 & 7 )))) v > . 
On the other hand as (u> 7 , ip~) is a fixed point of Q-y, we have that 

(>(™ 7 ,vg,(™ 7 ,^(i-G^ 7 )) 



(3.28) 



(3.29) 



Ij 7 ^,^ ) 



7/ up 



IKf + T^V (I-G'sWJ) . (3.30) 



2Wi 

It follows from (j2~g|) . and similarly to (|3~2"gj) . on noting (pOTj) and (j3~2l)]) that 



+ ^7 : ( J -^(^ 7 )) 



> 



4Wi 



[lhM" 2d ] 



> min 1 



4Wi /j, h -y 



9 ed\V\ 

y !— !■. 3.31 

' 2Wi V y 



Therefore on combining (|3.30j) and ()3.3ip . we have for all 7 sufficiently large that 

(>(™ 7 ,^ 7 ), (« 7J ^ (I- G' s (^)))) d < 0, (3.32) 

which obviously contradicts (|3.29[) . Hence the mapping T has a zero in B 1 for 7 
sufficiently large. □ 

Theorem 3.1. For any 8 G (0, \], > 1 and any partitioning of [0,T] into 
time steps, then there exists a solution {(Ug h , c^)}^ e \Yh x S"] Wt to (Pf{). 
In addition, it follows for n = 1, . . . , Nt that 



Re||u? fc - u^f + (1 - e)At m || Vix^JI 2 



2Wr 



2 Ai m / tr(/3 5 «J + [ft^)]" 1 - 21) 
m=l 

1 , p ™ 

^ ' m—l 

Moreover, it follows that 

KJI + ^IIKfcl-lO+E^ / H^KJ^II^c- 

(3.34) 



max 

n— 0,...,Nt Jt) 



n=l 



I 



I 
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Proof. Existence and the stability result Q3.33P follow immediately from Proposi- 
tions |U and ED] respectively, on noting pT20|) . (|3.11b|) . (|3.10a|) and ([2~T9]) . The 
bounds (pTBl)) follow immediately from (13"33")) . on noting (|2.7bl) . (|2~g|) . (fL7| and the 

fact that e R d S p d D for any e R d s xd . □ 



3.5. Convergence of (Pf* h ) to fP At J 

We now consider the corresponding direct finite element approximation of (P), i.e. 
(P A 1) without the regularization 5: 

(P At ) Given initial conditions (m£,<t£) e x S£ with er° satisfying (|3.11al b), 
then for n = 1, . . . , Nt find (it^ 1 , cr^) g V° x S° such that for any test functions 

K0)eV°xS° 



V L 



Re 



At r . 



Re 



"+t[(K _1 ^W)'"-<'(( 



■u 



n-l 



(1 - e)Vul : V« + — *l : V« 



At n 



0-2((V<X):0- 



Wi 



3=1 



kr 1 • n| Itrl: 







+n. 



= 0. 



We introduce also the unregularised free energy 



2Wi 



tr(</> - G(0) - I) , 



(3.35a) 



(3.35b) 



(3.36) 



v 



which is well defined for (v,<f>) E V° x S° with <p being positive definite on T>. 

Theorem 3.2. For all regular partitionings Th ofT> into simplices {Kk} k *i a nd all 
partitionings {At n }^ 1 o/[0,T], there exists a subsequence {{(it^ ^ CT 5 /J} ra =i}(5>ci> 
where {(u^, tr^)}^ E [V° x S°]^ SO foe S fP A <J, and {(ug/o-JJ)}^ € [V° x 
S^] such that for the subsequence 



u 



S.h 



5 ->• 0j 



/or n = 1, 



,N 7 



(3.37) 



/n addition, for n — 1, . . . , Nt, Ik^G Rjpj, fc = 1, . . . , Nk,- Moreover, 
{{<,<)} n=i G K x S^]^ sofoe S fPfj and for n = 1, . . . , N T 



F{ut,<Tl)-F{ul-\*l^) 



Re 



n-l 1 1 2 



At„ 2At n J v 

+ ( trK + K]- 1 -2/) 

2Wi Jt> 



(1-s) \\Vu 



D 



<- 2 (ls) 



ni|2 



1 + Cp 
2(1 -e) 



\f n \\ 2 H- 1 (V) 



(3.38) 



i 



i 
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Proof. For any integer n S [1,-ZVr], the desired subsequence convergence result 
(13.371) follows immediately from (I3.34p . as (ug h , <Tg h ) are finite dimensional for 
fixed V^xSg. It also follows from (f3734|) . (pH?) and ' ([217) that [<r£]_ vanishes on 
£>, so that er? must be non-negative definite on V. Hence on noting this, (|2.17j) and 
(|3.37|) . we have the following subsequence convergence results 

/3 5 (er£) ^ er£ as 5^0+ and As^/J a h as 5 ~> °+ • ( 3 - 39 ) 

It also follows from ([333]) . ([33^ and as [^(cr^)]" 1 /?^^™^) = J that the 
following subsequence result 

[ft (*?,,)] _1 -»• WhT 1 as ,5^0+ (3.40) 



holds, and so erJJ is positive definite on T>. Therefore, we have from (|3.37[) and (|2.ip 
that 

^K^GW) as ^0+. (3.41) 

Since w^ 1 , m^ 1 6 C*(P), it follows from the S£ version of IpTS]) . (l3~6l) and 
(f3737|) that for j = 1, ... , jV B and for all <j> e S° h 



N E 

E 



"L 1 • « 



-+-E/ K _1 -^K : 4> +url = E / K'^HKix-^^" 1 

k=1 JdK k j=l jE 3 

as <5 ^ 0+ . (3.42) 

Hence using (|3~37|) . ([3~3"5j) and (|3~4^|) . we can pass to the limit S -> 0+ in (P££), 
(l3~T2il b). to show that {«,<)}^i G [V° x S£] Wt solves (P^ 4 ), (l3~35a b). Sim- 
ilarly, using (j3~37| . (j3~39)l . (jOU)) and (pT4T]l . and noting (|2~20|) and (|3~36)) . we can 
pass to the limit S — > 0+ in (|3 . 15[) to obtain the desired result (I3.38[) . □ 

Remark 3.1. Most numerical approximations of (P) suffer from instabilities when 
Wi is relatively large, the so-called high Weissenberg number problem (HWNP). 
This problem is still not fully understo od. Some reasons for these instabilities 
are discussed in Boyaval et a/.j I e.g. poor numerical scheme or the la ck 

of existence of a solution to (P) itself. In addition in Boyaval et a ;J BLM0 9l g_ 
nite element approximations of (P) such as (P^*), approximating the primitive 
variables (u,p, <x), are compared with finite element approximations of the log- 
formulation of (P), introduced in Fattal and KupfermanJI-ESI which is based on 
the variables (u,p, ip), where ip — In <x. The equivalent free energy estimate for this 
log-formulation is based on testing the Navier-Stokes equation with u as before, 
but the log- form of the stress equation with (exp-0 — 7). Whereas the free energy 
estimate for (P) requires <r to be positive definite, due to the testing with lner, the 
free ener gy estima te for the log- formulation requires no such constraint. In Boy- 
aval et all BLM09 | a constraint, based on the initial data, was required on the time 
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step in order to ensure that the approximation to a remained positive definite for 
schemes such as (P^ *) approximating (P); whereas existence of a solution to finite 
element approximations of the log-formulation, and satisfying a discrete log-form of 
the free energy estimate , were shown for any choice of time step. It was suggested 
in Boyaval et a /| BLM0 9] ^.j^ ma y k e ^he reason why the approximations of the 
log- formulation are reported to be more stable than those based on (P). However, 
Theorem 13.21 above shows that there does exist (at least) one solution to (P^*), 
which satisfies the free energy estimate (|3.38[) . whatever the time step. Of course, 
we do not have a uniqueness proof for (P^*). 



4. Regularized problems with stress diffusion and possibly the 
cut-off (3 L 

4.1. Regularizations , (P$P' ), of (P) with stress diffusion and 
possibly the cut-off f3 

In this section, we consider the following modified versions of (P) for given constants 
a G JR >0 and L > 2: 

(Pi L) )Find «i L) : {t,x) G %T)xV^v!£\t,x) G M. d , p^ : (t,x) G (0,T)x2?^ 
pi L \t,x) G K and o- ( a L) : (t,x) G [0,T) x V i-4 a^\t,x) G R d s xd such that 

Re + («W ■ V)«^ = - VpW + (1 - e)Au^ + JL diy^^W) 

+ / on £> T , (4.1a) 
div«W=Q onV T , (4.1b) 



tr^ aA<ri L) on V T , (4.1c) 



dt 

1 

~ Wi 

u { a L) (0,x) = u°(x) VxeV, (4.1d) 

o-W (0, x) = cr°(x) VxeV, (4.1c) 

u^=0 on(0,T)xd£>, (4. If) 

(n ro • V) CT W =0 on (0, T) x dV ■ (4.1g) 

where nx> is normal to the boundary dT>. 

Hence problem (P^) is the same as (P), but with the added diffusion term 
aAff^' for the stress equation (|4. lc|) . and the associated Neumann boundary con- 
dition ( |4.1gP ; and in the case of (P^) with certain terms in (|4.1a[ c) involving er^ 
replaced by /3 L (<x^), recall (|2.3[) . Of course, it is naturally assumed in (P^) that 
er^ is positive definite on T>t in order for f3( L \cr^) to be well defined. 

We will also be interested in the corresponding regularization (P^ L j) of (P^) 
with solution (u^\,p^ s , f^]); where ^ L \-) in (|4.1ab g) is replaced by (3$ (■), and 
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so that <x^ <s is not required to be positive definite. 



4.2. Formal energy estimates for (P^l) 



Let F^iu^j.a^l) denote the free energy of the solution (w^, P^ a )a,s) to 
problem (P^$), where Fg : W x S — >• M is defined as 



(L) JL) ^(i). 



Re 



/jWI 2 + ^X^-Gf>W-/) 



2 J v " " 2Wi 
We have the following analogue of Proposition 12. II 



(4.2) 



Proposition 4.1. Let , , o~^\) be a sufficiently smooth solution to problem 



Q?a})- Then the f ree energy ' (u^j , er^j) satisfies for a.a. t G (0,T) 



d 
It 



—F {L) (u {L) <t {l) ) 



IVu 



( L )||2 



+ 



2Wi 2 



tr(/3i L) (^]) + [/3i L) (TS)]- 1 -2/) 



aecr 
2WT 



(4.3) 



Proof. Multiplying the Navier-Stokes equation (|4.1al) with and the stress 

equation ()4.1c[) with ^[(I ~ G^ (c^])), summing and integrating over V yields, 
after using integrations by parts and the incompressibility property in the standard 
way, that 



v 



Red_ 
~2~di [ 



u 



< £ >" a + (l-e)||V«2|| a 



£ 

Wi Jv i 
1 ffd 



PS \ a a^) ■ Vt V<5 



2 vv ^"a,^«i-v)^r«i) 
^((v w S)# ) ( CT S) + 4 L) (-S)(v<i) 



Wi V 

T 



Similarly to (|2.12[) . we have that 



Ver 



(L) 



a.e. in 2?y. 



Using (12.16[) . we have that 



-<T {L) 

dt a ' s 



(4.4) 
(4.5) 
(4.6) 
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We will deal with the convection term differently to the approach used in (j2.23[) . 
as that cannot be mimicked at a discrete level using continuous piecewise linear 
elements to approximate c^J. Note that we cannot use S® with the desirable 

property (|3.5[) to approximate c^ L j, as we now have the added diffusion term. 
Instead, as er^j has been replaced by (it^]) = (G^ L ' (c^])), on recalling 
(|2.6|) . in this convective term and as u^l G V, we have that 



= //3f ) (-S):(«S-V)Gf)V^ 
Jt> 



= jy"l-V)tv(Hf\Gf ) \ ( T^)))=Q 1 (4.7) 



where we have noted the spatial counterpart of (|2.16[) . Similarly to (|2.24[) and (|2.25l) 
we obtain that 

= 2tr((Vu^^ (4.8) 



and once again the terms involving the left-hand side of (|4.8[) in (|4.4|) cancel with 



the term (o^j) : Vu^] in (|4.4[) arising from the Navier-Stokes equation. 



Finally, the treatment of the remaining term (tr^g — I^j '■ (i — G S L ** (c^j)) follows 
similarly to (|2.26[) : and so we obtain the desired free energy inequality (|4.3p . □ 



The following Corollary follows from (|4.3|) on noting the proof of Corollary [2TTJ 

Corollary 4.1. Let {u^\,p^\,<J^\) be a sufficiently smooth solution to problem 
(P^ L ]). Then it follows that 

2 



sup Ff)(«g(t I .),«rS(t,.)) + ^ / llVGfVS)!! 2 
te(o,T) ^ Wl JD T 



(1 - £ )||V<il! 2 + r^tr(/3 2) + t/3r J (< )]^ - ^) 
Wi 

2(1 -e) 



< 2 ( *r w°) + ^tt^s \\ni H0 , T , H -, m i . (4.9) 



5. Finite element approximation of (P„g) and (P^) 
5.1. Finite element discretization 

We now introduce a conforming finite element discretization of (P^J), which sat- 
isfies a discrete analogue of (|4.3[) . As noted in the proof of Proposition @TT] above, 
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we cannot use 5° with the desirable property (|3.5[) to approximate er^J , as we now 
have the added diffusion term. In the following, we choose 

:= W 2 h c W or W^' + c W , (5.1a) 

Q\ = {qeC{V) : q\ Kk e¥ 1 k = l,...,N K }cQ, (5.1b) 

S{ = {0 G [C(V)]f d ■■ <t> \i< k e [Pijf d k=l,...,N K }cS (5.1c) 

and Vl = (v£Wl: J gdivi> = Vg G Q{\ ; (5. Id) 

where is defined as in (|3.3a|) and, on recalling the barycentric coordinate nota- 
tion used in (13.41). 



W*'+ := «e[C(D)fnw : v \ Kk e 



ispanJJ?7. 



k 

>l 

i=0 



k = l,...,N. 



K 



(5.2) 

The velocity-pressure choice, W? x Qj, is the lowest order Taylor-Hood element. 
It satisfies (I3.1[) with W° and Q° replaced by and Q^, respectively, provided, 
in addition to {Th}h>o being a regular family of mes hes, tha t each simplex has at 
leas t one v ertex in X>, see pl77 in Girault and 

RaviartPBS] 

in the case d = 2 and 

B gjBoffl7j j n cage ^ _ 3 Qf coursej this is a very mild restriction on {Th}h>o- 

The velocity-pressure choice, W^' + x Q^, is called the mini-element. It satisfies (|3.1[) 

with W° and Q° replaced by W! ' + and Q^, respectively; see Chapter II, Section 

4.1 in Gir ault an d Raviari! ^ in the case d = 2 and Section 4.2.4 in Ern and 
GuermondESQl 

in the case d — 3. Hence for both choices of W\, it follows that 
for all v G V there exists a sequence {vh}h>o, with G V^, such that 

lim \\v - v h \\ H im) =0. (5.3) 

We recall the well-known local inverse inequality for 

\\q\\ La o (Kk) <C\K k \- X [ \q\ VqeQl k = 1, . . . , N K 
JK k 

=> \\x\\L~(K h ) ^C'lK,]- 1 f llxH V X eSl k = l,...,N K . (5.4) 

We recall a similar well-known local inverse inequality for 

\\Vv\\ LHKk) <Ch^\\v\\ L 2 {Kk) VveVl k = l,...,N K . (5.5) 

We introduce the interpolation operator 717, : C(T>) — > Q^, and extended natu- 
rally to TT h : [C(V)]f- d -> Sj,, such that for all rj G C(P) and G [C(©)]f xd 

Tfc^C-Pp) = V{P P ) and TTh^Pp) = 0(P P ) p = l,...,iVp, (5.6) 

where {P P }p=i are the vertices of Th- As 6 S^ does not imply that G^ L ' (0) G S^, 
we have to test the finite element approximation of the (P„j) version of (|4.1cj) 
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with I — nh[Gg' (cTaf^h)] G S^, where tr^g^ £ is our finite element approx- 
imation to tr^l at time level t n . This approximation of the (P„ J) version of 
(|4.1c[) has to be constructed to mimic the results (|4.5|) - (|4.8|) . when tested with 

/-^[Gj L) V£S)]eSi. 

In order to mimic (|4.5|) we shall assume from now on that the family of meshes, 
{Th}h>o, for the polytope T> consists of non-obtuse simplices only, i.e. all dihedral 
angles of any simplex in Tk are less than or equal to \. Of course, the construction 
of such a non-obtuse mesh in the case d = 3 is not straightforward for a general 
polytope V. We then have the following result. 

Lemma 5.1. Let g G C 0,1 (R) be monotonically increasing with Lipschitz constant 
<7Lip- As Th consists of only non-obtuse simplices, then we have for all q G Q^, 
<p G that 

9u P VTT h [g(q)] ■ Vq > \\Vn h [g(q)}\\ 2 and g Lip Vv h [g (</>)] :: V0 > || V7r h [g(cj))}\\ 2 

onK k , k = l,...,N K . (5.7) 

Proof. Let K\. have vertices {Pj }f = o, and let r)j(x) be the basis functions on Kk 
associated with and S^, i.e. rjj |^ fc e Pi and rn{P^) = Sij, i,j = 0, . . . ,d. As 
Kk is non-obtuse it follows that 

Vryf • Vrf < on K k , i,j = 0,...,d, with i^j. (5.8) 

We note that 

d 

E Vj = 1 on K k => 

j=o 

d 

\\VVif = - E • V//;' on K k , i = 0,...,d. (5.9) 

3=0, 



Hence for a,, 6,- G R, i = 0, . . . d, we have that 

v(E^)-v(]Tv/*) = E 



i=0 3=0 i=0 



d 

a<MV>7?H a + E "' /, - v '/; • V//^ 



= -E E ai^i-^v^.v^* 

= - E E( a * - ( 6 « - h i) v ^ ■ Vr '" • ( 5 - 10 ) 

i—0 j>i 

Similarly for a;, bi G M^ xd , i = 0, . . . , d, we have that 

d d d d 

V(E <* ^) :: V(E ^ r?|) = - EE " ^ : ( b > " M V ^ ' V? >* ■ ( 5 ' U ) 

i—0 j=0 i=0 3>i 
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The desired result (|5.7[) then follows immediately from (|5.10|) . (|5.11j) . (|5.8j) and our 
assumptions on g. □ 

In order to mimic (|4.6[) and (|4.8[) . we need to use numerical integration (vertex 
sampling). We note the following results. As the basis functions associated with 
and Si are nonnegative and sum to unity everywhere, we have for k = 1, . . . , Nk 
that 

< ^ h [U\\ 2 ]){x) VxeK k , V</>G [C(T^)] dxd . (5.12) 
In addition, we have for k = 1, . . . , Nk that 

/ M 2 < / 7r h [\\x\\ 2 ]<C [ Hxll 2 VxeSj,. (5.13) 

The first inequality in (|5.13[) follows immediately from (I5.12[) . and the second from 
applying (|5.4j) and a Cauchy-Schwarz inequality. 

In order to mimic (|4.7[) . we have to carefully construct our hnite element ap- 
proximation of the convective term in the (P^fj) version of (|4.1cp Our construction 
is a non-trivial extension of an approach that has been used in the finite element 
approximation of fourth-order degenerate n onlinea r parabolic equations, such as the 
thin film equation; see e.g. Griin and 

RumpPEQO] 

and Barrett and Numbers! I 
Let {ei}f =1 be the orthonormal vectors in M. d , such that the j th component of 
is Sij, i,j = 1, ...,d. Let K be the standard open reference simplex in R d with 
vertices {Pi}f =0 , where Pq is the origin and p = e,-, i = 1, . . . , d. Given a simplex 
Kk G Th with vertices {Pj }f = cu then there exists a non-singular matrix B k such 
that the linear mapping 

B k :x eR d ^ Po + B k x G R d (5.14) 

maps vertex Pi to vertex P*, i = 0, . . . , d. Hence Bk maps K to Kk- For all r\ G 
and Kk <ETh, we define 

rj{x) :=r)(B k (x)) Mx g K => X?r)(B k (x)) = {BD^V^x) VS g K , 

(5.15) 

where for &\\xe K 

[^v(x)] j = ^v(x)=v(P j )-v(Po)=v(Pj)-v(Po) j = !,...,£ (5-16) 
Such notation is easily extended to </> G 



Given </> G and K k G %, then first, for j = 1, . . . , d, we find A^ L) ((/>) G Ms Xd , 



which depends continuously on <f>, such that 

: 4^ [ ^ L) ' ( ^ )] = s|:^ [tr( ^ i)(G * L 3 ' ($)))1 ° n (5 ' 17) 

where (nhrf)(x) = (irh7])(Bkx) for all a; G K and 77 G C(Kk). This leads to a unique 
igj (4>) . For the construction of 



choice of A^i- (0). For the construction of A.gh<f>) in the simpler scalar case (d = 1), 



1 
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see p329 in Barrett and NiirnbergPSQH To construct A s L j((p) satisfying (|5.17p . we 
note the following. We have from (|2"3|) . and (|2TT5]) that 

(3 { 5 L \cf>(P*)) : (Gi L) '(0(^)) -G< L) '(0(P o fe ))) 

<tr{Hf\Gf ) \<l>{P*))-Hf\Gf ) \ ( t>{P*))) 

< PPWP*)) : (G^'(0(Pf)) - G™' (ft!*))). (5.18) 

Next we note from ([23]) . (f2TB]l . ([2~7f) and ([i~2b| that 

- - /3 A L (0(P o fe ))) : (<#' (0(Pf)) - G L S ' W>(P fe ))) 

> i- 2 (0(Pf )) - /3 5 L (0(P o fc ))|| 2 ; (5.19) 

and so the left-hand side is zero if and only if fig (<p{P^)) = (3$ (4>(Pq)) . Similarly, 
we see from (|2"3]) . (EU) and the proof of (pjf) : that is, (j2~T3"|) ; that 

- (ft(0(P/)) - (3 S (4>(P k ))) : (G^Pj 5 )) - G' 5 (0(P o fc ))) > (5.20) 

with equality if and only if [3 s {(j>{P k )) = /3 s (<p(P$)). Hence, on noting ([536)1 . ([B38| . 
(|5TT9|) . (|530|) and (|i~2b)l . we have that 

lf]Q>) := (1 - Xf])pf\^)) + \f]pf\^)) 

if (^^(P/)) - fif \cj>(P k ))) : (Gi L) '(0(P/)) - G^Wo ))) ^ , 

(5.21a) 

A$(£) == pf ) (0(^i )) = /3f > Witf)) 

if (/3i L) (0(P/)) - /?f } (0(P o fc ))) : (G< L) '(0(^)) - G< L) Wo ))) = 

(5.21b) 

satisfies (|5.17l) for j = 1, . . . , d; where Xgj G [0, 1] is defined as 
'ttiHfHGW'MP? )) Pi L) (G< L) '(0(P o fc ))) 

-Pf } (0(Pf)) : (G< L) '(0(Pf)) -Gf >W<?)))' 



A 5j — 



(/3< L) (<KP fc )) - ^W/))) : (Gf } (0(Pf )) - Gi L) (0(P o fc ))) 

Furthermore, Ag L j(4>) € Mg Xd , j = 1, . . . ,d, depends continuously on <p |kv 
Therefore given G S^, we introduce, for m,p — 1, . . . , d, 

d 

A £, P (<£) = EK fl * [P, T ] JP e Mf d on *r fc , 

3=1 



fc = l,...,JV x . (5.22) 



It follows from (|5~2"2"|) . (j5~TT)) and ([535]) that 



A& P (0)~/3l L) W>)<W m,p=l,...,d; (5.23) 
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and for m — 1 , . . . , d 

X>2» : -Sr^^'iW = -£-*hMh¥\g£ } '(0)))] on K k , 



p=l 



' p 



k = l,...,N K . (5.24) 

For a more precise version of (|5.23l) . see Lemma [5.31 below. Finally, as the parti- 
tioning 7~h consists of regular simplices, we have that 

KBlr'WWBlWKC, k = l,...,N K . (5.25) 

Hence, it follows from ([Q2]) , ([535)) . (|5.21al b) and (|2~4| that 

l|Af; m , p (0)|U»( D ) <CL V0eS,\. (5.26) 

5.2. j4 /ree energy preserving approximation, (Pa'sh*)' °f (^ >< al) 

In addition to the assumptions on the finite element discretization stated in sub- 
section [5TTJ and our definition of At in subsection 13.11 we shall assume that there 
exists a C e M>o such that 

Ai„ < CAi„_i, n = 2,...,N, as Ai ^ 0+. (5.27) 

With Ati and C as above, let Ato G M>o be such that Ati < CAt$. Given initial 
data satisfying (|2.19[) . we choose it° G and cr° G throughout the rest of this 
paper such that 

/ [u° h ■ v + Ai Vu£ : Vv] = [ u° ■ v Vu G , (5.28a) 

/ : X ] + AtoVrrO :: V X ] = f <r° : X V X G S£ . (5.28b) 

It follows from (j5.28al b). (|5TT5j) and (|2TTT?1) that 

/ [Kl| 2 + Kl| 2 + At [||V<|| 2 + ||V<t°|| 2 ]] <C. (5.29) 
In addition, we note the following result. 
Lemma 5.2. For p = 1, . . . , TVp we ftaue i/iai 

^LJ£|| 2 <^(^K<^cII£|| 2 V^eR d . (5.30) 
Proof. It follows from (|5.28bj) that 

/ [TTfcKo-g - <x° lin 7) : x] + A* V(<r° - a^J) :: V X ] = / (<r° - <r° min /) : * 

V X eS£. (5.31) 

Choosing x = $,$, T i], with f? G yields that := £ T (ff° — <r^ in J)£ G satisfies 
/ [n h [z h 7]} + At Vz h -Vr]]= [ z V V?/ e Qf , (5.32) 
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where z := £ T (<x° — cr^^J)^ G L°°(T>) and is non- negative on recalling (|2.19p . 

Choosing r\ = 7r^[z/i]_ G Q^, it follows, on noting the version of (|5.13p and 
(|577|) with g(-) = [•]_, that 

/ [^[^-P + AioHVTT^]-!! 2 ] < / [7T h [[Zft] 2 .] + At Vz h - V-K h [z h }-] 
Jv Jv 

= [ zn h [z h ]- <0. (5.33) 
Jv 

Hence nh[zh]- = and so the first inequality in (|5.30p holds. Repeating the above 
with er^ijj and [•]_ replaced by <x° lax and [•] + , respectively, yields the second 
inequality in (|5.30p . □ 

Our approximation (Pl^^ ) °^ (Pal) * s then: 

(PSD Setting (<j;>SJ) = «, <r°) G Y\ x S*, then for n = 1, . . . , N T 
find (tifjj, c^sh) £ Vj, x Sj such that for any test functions (v, <f>) G x 



<L,)n-l> 



Re ( 

Re 
~2~ 



Atr, 



(.f»^(2>), (5.34a) 



,5,h 



r (L,)n-l- 
a.S.h 



At r 



^ m— 1 p— 1 P 

In deriving (P^f *), we have noted (|3TT5]) , (H~4al) and (j5~2^)) . 



(5.34b) 



Before proving existence of a solution to (P^/f* ), we first derive a discrete 
analogue of the energy estimate (|4.3p for (P^fj). 



5.3. Energy estimate 

On setting 



Re 



V 



S 

2Wi ./.p 



J^ h [tr(4>-Gf\4>)-l)] 



(5.35) 



we have the following discrete analogue of Proposition 14.11 
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Proposition 5.1. For n = 1,...,Nt, a solution (i^gh^a'sh^j e x ^° 
{5.34a\ b), if it exists, satisfies 

t?{L) f (L.)n (L,)n, „(£) , (L,)n-1 (L,)n-U 

^S,h \ U a,S,h i a a,S,h ) ^ S,h \ U a,5,h ^a.S.h ) , He / (£,)„ (i,)n-l,,2 



+ (1 -e) £ ||V<^;;|| 2 + ^^„[tr( / 3< i )(<755) + [^'(.tSD]- 1 - 2/)] 



< ^1- s) jf ||V^>f + ±0L\\rf H - Hv y (5.36) 



Proof. The proof is similar to that of Proposition 13. 1[ we choose as test functions 
« = 6 VJ and = ^ (l-TTfcl^'C^)]) 6 S£ in »), and 

obtain, on noting (|3.14p . (|2.7al d.e). (|5.7p with g = —Gg having Lipschitz constant 
S- 2 , dnHU» and (|Q5| that 
/ f « „ (-k.) n \ 



Tp{L) ( (L,)n (L,)ns_ „(L) , (L,)n-\ (L,)n-U 

>- ■— ^ ' + (l-e)/ ||V« 



(i,)™ ||2 

a,S,h II 

2? 



Re f „ (£,)„ (i,)n-l||2 



2A*„ 7 P " aAh a ' s ' h " ' 2Wi Jr 



2Wi 2 



+ / -SIT 1 • VTT^tr^f J'(.r£>)))] • (5.37) 

The first desired inequality in (|5.36|) follows immediately from (|5.37|) on noting 
(|5Td|l . (|3~2a| . (|L9j) and that n h : (7(2?) -> Q^. The second inequality in (fOg]) 
follows immediately from (12.281) with i/ 2 = (1 - e)/(l + Cp). □ 

5.4. Existence of discrete solutions 

Proposition 5.2. Given (u^g? l i (T ^a5h *) e ^ x ^ft and for any time step 
At n > 0, then there exists at least one solution ^ti^, ^aSh^j ^ V\ x ^° 

low . 



Proof. The proof is similar to that of Proposition 13.21 We introduce the following 
inner product on the Hilbert space x 

((to, </>), (v, 4>)) h v = [ [wv + n h [xp : 0]] V(u>, </>), (v, 0) G V,\ x S x h . 
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Given r« (£ ' )n_1 o-^')™- 1 ^ 
uiven ^« ai5) /, , <J a ,6,h 

any tfi) e V[ x S,\ 



G Vi x S^, let P l : Vl x Si -> Vi x Si be such that for 



(^(to,^),^,^)) 

r /« 

Re - 



(L,)n— 1 ' 
a,8.h 



• « + (1 - e)Vic : Vv + ±-Tr h [^ s L \xp)] : V« 



Re r 
~2 



(< 



(L,)n—1 



aVip :: V</> 



2Vti;: 7 r /l [0/3f ) (^)]+7r /l 



d d 



^ m— 1 p=l P 

V(«»eVixSj. (5.38) 

A solution (u^a8h' a< a5h) ^° Q3.12a[ b). if it exists, corresponds to a zero of J 7 ' 1 . 
On recalling (|5.22[) and (|5.21a[ b). it is easily deduced that the mapping T h is 
continuous. 

For any (w, ip) G V£ x S^, on choosing {v, cj>) = (w, [I - k h [G ( s L) (t/>) 
we obtain analogously to (|5.36|) that 



(£), (L,)n-l (L,)n-U 



> 



1 -£ 

2 

aeS 2 



At n 



2Wi 2 



Re 
2Ai n j-j? 



,(i,)"-l||2 
1-1 



I + Cp 



^ jf ||V,r h [Gf >'(*)]|| a - ||r|| 2 ff _ 1(P) . (5.39) 



Let 



(v,4>)\\ h D := [((«,0),(«,0)) 



P 



[IMI 2 + **[IMI 3 ]] 



If for any 7 G R>o, the continuous mapping J- h has no zero (w^r CT^Pl 1 ), which 
lies in the ball 

B*:={(t7,0)GVixSjt : ||K0)||^< 7 } 5 

then for such 7, we can define the continuous mapping Q^-.B^^t B^ such that for 
all (v, (j>) G 

J"' l (v,0) 



l^(«,0)ll 



I 



I 
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By the Brouwer fixed point theorem, C7 7 has at least one fixed point (iu 7 ,7/> 7 ) in 



S 7 . Hence it satisfies 



||(^ 7 ,^)|£=||^K,V 7 )||i=7- (5-40) 
On noting (15. 4j) . we have that there exists a /i/, € M>o such that for all G Si, 

(5.41) 



li = (X>) < lkft[||0ll <Vh j^h[\\4>\\ 

It follows from (|535|) . (|2~8j) . (pTITJ) and (f^40|) that 



> 



> 



7 

Re 
~2~ 
Re 

T 

Re 



2Wi 



4Wi 



7r h [tr(^ 7 -G^(^)-I)] 



4Wi/^7 



.'Re 
> mm — 



2 '4Wi^ 7 y \J V 
Re e 



/ 7T fc [||^ 7 ||]-2d|2?| 

R[||V 7 ll]IU~p>) 

r,, -^|| 2 +7r h [||^J| a 



gd|P| 
" 2Wi 



2 '4Wi/x h7 



7 



ed\V\ 
2Wi 



Hence for all 7 sufficiently large, it follows from (|5.39|) and (|5.42j) that 



(j-' l ( U ; 7 ,V 7 ),(^ 7 ,^(/- 7 r,[Gi L) '^ 7 )]))) 



> 0. 



v 



ed\V\ 
2Wi 



(5.42) 



(5.43) 



On the other hand as (u> 7 , i/O is a fixed point of C/' 1 , we have that 



(^ l ( U ; 7 ,^ 7 ),( U , 7 ,^(/-^[4 L) '(^ 7 )]))) 



I J r ' l (u; 7 , ip 



7' Y^JW-p 



IKH* + ^^[^7 : ( J - G i L) '^7))]] ■ ( 5 - 44 ) 



It follows from (f2T8]> . and similarly to f|5.42[) . on noting f|5.41[t and (|5.40[) that 



2Wi 



r7r fc [^: f-G™ WO ] 



> 



4Wi 



> min 1 



4Wi ^7 



Ml!<MN-2d] 

' 2Wi ' 



(5.45) 



Therefore on combining (|5.44j) and (15.451) . we have for all 7 sufficiently large that 



2Wi 



/-7r,[Gf ) '(V 7 )])))^<0, 



(5.46) 
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which obviously contradicts (|5.43l) . Hence the mapping T h has a zero in B 1 ^ for 7 
sufficiently large. □ 

We now have the analogue of stability Theorem 13. II 

Theorem 5.1. For any 5 G (0, 4], L > 2, Nt > 1 and any partitioning of [0,T] 

into Nt time steps, there exists a solution {(u^shi "'ijl)}^ ^ [Yh x ^/J"^ ^° 
/p(i,)At , 

In addition, it follows for n = 1, . . . , Nt that 



•Y(i)V_(i,)"»n|i2 
j,«,h Jill 



m=l 



m— 1 



< F$(ul<r° h ) + £ A^IIZ-H^.^ <C. (5.47) 

^ ' m=l 

Moreover, it follows that 

max / [ll^tl^ + ^IIkKin+^^tll^SS-ll] 

n=0,...,N T Jx> L - 

+ £ / AMVugJU 8 + At„7r h [ HDJW^)]- 1 !! ] + ||u£# - nSr'll 2 

<C. (5.48) 



Proof. Existence and the stability result (|5.47[) follow immediately from Proposi- 
tions E2 and EU respectively, on noting (f5T35|) . (f5~29]) . ([530]) . (|3.10a|) and (|2~T9]) . 
The bounds (f5T48|) follow immediately from (|5T47j> . on noting (|2.7bj) . (|2~8]) . (fl~7| 

and the fact that /^ L) (</>) E R d s x P d D for any G M^ xd . □ 

5.5. Convergence of (P%>£ l ) to (P^ ) 

We now consider the corresponding direct finite element approximation of (P^), 
i.e. (P^^ At ) without the regularization 5: 
We introduce 

Sh,PD = {<t> 6 S£ : 4>{P P ) G Rf P d D for p = l,...,N P }c S PD . (5.49) 
It follows from ([530]) that <x° G S^ P23 . 

(P^ )At ) Setting («g'VS°) = «,<t°) 6 Vi x S^, then for „ = 
1, . . . , Nt find («^'", ^ah*) £ x Sj, such that for any test functions (v, <p) G 



1 



1 
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Re 



u 



(L,)n 
cx.h 



U 



(L,)n-1 * 
a . h 



At r 



Re 
T 



((«£^" l -V)« 



(1 - £ )V<> : V* + ^P< £ > : V « 



= </ B ,«>Hg(i>), (5.50a) 



(L,)ra-1* 



+ « / " V*- 2 / V«£> : <r k [*^>(<>)] 

- / E £>ST ^ AW : ^ = . (5.50b) 

^ m— 1 p— 1 ^ 

Remark 5.1. Due to the presence of f}( L > in (|5.50a[ b). it is implicitly assumed that 
°"S" G S Ipd> n = l,...,iV T ; recall (& In addition, (0) for € PD 
is defined similarly to (|5 . 22f) with Ag L j(<p) replaced by A^ L \<p), which is defined 
similarly to l|53Ial b) with \$,$ L) and Gf replaced by \f\ (3^ and G( L \ 
with \f ] defined similarly to with pf\ Gf ) and i^ L) replaced by ^ L \G^ 
and H( L \ Hence, similarly to (|5 .26[) . we have that 

||A^ p (0)|Uoo (P) <Ci V0GSi, PD . 

Wc introduce also the unregularised free energy 



v,<j>) 



Re 

T 



2Wi 



7T/, 



[tv^-G^ (</>)- 1)}, 



(5.51) 
(5.52) 



which is well defined for all (v,<p) £ V\ x PD . 



Theorem 5.2. For aZZ regular partitioning s Th of T> into simplices {Kk}u=i and 
all partitionings {Af n } n ^j o/ [0, T], there exists a sub- 
seauenee {{(«£> ^^^o, ufere {(ug£, crg*)}^ E [V* X Sj,]^ 

sofoes ^Sf'A {( M S'>£ )n )}n=i G [Vfc x Sfe]^ stlc/l ^ai/or tte sufese- 
owence 



:,5,h 



r (i.)r 

!,/» 



In addition, for n = 1, 



J /i.P_D 



as <5 -> 0+ , for n=l,...,N T . (5.53) 

h,P.D) 



.*r, £ Sl PD , and {(«<$V&*)fe £ [Y£ X 



Moreover, we have for n = 1, . . . , iVy £/iai 

L (L,), 



p{L),(L.)n (L,)n, „(L) , (L,)n-1 (L,)n-U 
b h \ U a,h > CT a,h ( M a.ft > CT a.h ) 



At* 



Re 
2AT 



n JT> 



,(£,)'■ 



(L,)n-1,,2 



I 



I 
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+ (1 - e) [ \\Vu™»\\' + ^~s[ ^[tr(/?V^ )n ) + ^{fr^TT 1 2/)] 
Jv 2Wi Jv 



IV ^VVl Jv 

<\{l-e) jjVu^f + i±^L||r||? r - 1(1 , )> (5.54) 



and 



max 

n— 0,...,Nt Jd 



(L,)n<\2 , r II _(-£-. ) n l 



• £ / [At„||v<f i^ + At^iii^^g")]- 1 !!; 



(L,)n (i,)n-l||2 



<C. (5.55) 



Proof. For any integer n 6 [1, iVr], the desired subsequence convergence result 
f|5 .53[) follows immediately from (|5.48l) . as (u^')2, ^a Sh) are nm te dimensional 
for fixed VjxSj. It also follows from (JEHJ), (15331) and (l2~T7l) that 7r h [ [er^]_] 
vanishes on D, so that must be non-negative definite on 2?. Hence on noting 

this, (|2.3[) . (|2 . 1 T[) and (|5 -53[) . we have the following subsequence convergence results 

/fViT)^^") and VS?D as ^° + - 

(5.56) 

It also follows from (ET481) . (l536| and as [^VSDI^^KD = J that 
the following subsequence result 

holds, and so g S,\ P£) . Therefore, we have from (15.531) and (|2.1I) that 

Gt\^ h )^G^{a^ n ) as ^ 0+ . (5.58) 

Similarly to (15351) . it follows from (15331 . (15361 . (15T221 and <(53I5I b) as er^'™ e 
P£3 that for m, p = 1, . . . , d 

aS>55) A ^(^S n ) as « -><»+■ ( 5 - 59 ) 

Hence using (I5.53[) . (|5.56[) and (15.59[) . we can pass to the limit 5 — >• + in (P^^f*), 
»), to show that {(uS n ,^>)}^ £ [Vi x S^]^ solves (P^ A <), 
(|5.50a[ b). Similarly, using ([5T53]) . ([53B1) . ([537]) and ([B35]) . and noting (fOS]) and 
()5.52p . we can pass to the limit S — > 0+ in (|5.36j) and (|5.48|) to obtain the desired 
results (|5.54|) and (|5.55|) . □ 



For later purposes, we introduce the following notation in line with (|3.9[) . Let 

l( ah At e cd ' T l; Y i) and M i L ^ )At,± e T ; v h) be such that for n = i, ... , Nt 



<t% ') := L A^S )n (0 + ^"SMO * e (5.60a) 



I 



I 
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(L,)At,+ 
l a,h 



(V) :=«&?"(•) 



,(i,)Ai, 



. a>fc ■(t,0:=«S"" 1 (0 (5.60b) 

(efr 1 ,!"), (5.60c) 



and A(t) := Ai„ 
We note that 



(L,) At 



(L,)At (i,)At,±_ r . 



0t 



te(i B-1 ,t n ), n = l,...,JV T , (5.61) 



where t\ := f" and t n _ := t"- 1 . We define <r^ )A * G C([0, T]; Sj l PD ) and 

°S? At,± e L°°{Q,T- Slpjj) similarly to jUg b). 

Using the notation (|5.60al b). (|5.50al) multiplied by At n and summed for n = 
1, . . . , iVr can be restated as: 

■ Bn (L ' )At 1 
Re^l^ •» + (!- e)V«i L , )A ^ + : V« <ft 



o ./P 



Re 



2 Jo 7x> 
f T r 



(L,)At,— — , (L.)At.+ , (L.) At.— 



(L,)At,+ 



rff 



Vv G L 2 (0,T;V{). (5.62) 

Similarly, (I5.50bp multiplied by At n and summed for n = 1, . . . , Nt can be restated 



dt 



as: 



o Jx> 



- a 



5<T (i,)At 



0t 



Ver 



(L,)At,+ 



: Vxtft-2 



o Jx> 



(L,)At,+ 



J» 



^[x/? (L) (-i L , )At ' + )]^ 



„T „ d d 
'o JV 



E B<* 

m— 1 p=l 



(i,)At,-| a (L) ( (L,)At, + ^ . 



= 



V X ei 2 (0,T;Si). (5.63) 
We note also the following Lemma for later purposes. 
Lemma 5.3. For all Kk G 7h> an d f or a ^ G Sj PD we have that 



||^W(0)]-/3W(0)f + max 
Kfc m,p=l,...,rf 



Jmp | 



||AW(0)-^)(0)«5„ 

< Ch 2 f ||V0|| 2 . (5.64) 



Proof. First, we have from (f2TT7|) that for all 4> G 



PD 



3=0 
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< 



2 



C\K h \ ^||0(P/)- V |iL» ( x fe ) 



3=0 



< Ch 2 \K k \ ||V0||io. (Jri>) < Ch 2 f ||V0|| 2 . 

(5.65) 



where {Pj} d =0 are the vertices of K k . Hence we have the desired first bound in 
(15^41 . 

It follows from the 5 independent versions of (|5.22j) and (|5.21a[ b). recall Remark 
[S~T1 ([OB]) and ([2~T7| that for all 6 S£ PD 



l|2 

rnp 



||AW(0)-7T h [^)(0)]* 

ii £ [k b * [Af - MP (L) m ] [bi 

3=1 

<C [ J2\\Af^)-n h ^ L H4>)}\\ 2 



< C\K k \ . max ||/^>(0(P/)) - /3^(^(Pi))f 

1,3=0,. ...d J 

<C\K k \, max ||0(P/)-0(if)|| 2 

< C/i 2 / ||V0|| 2 . (5.66) 
Combining (|5.66|) and the first bound in (|5.64l) yields the second bound in (15.64|) D 

6. Convergence of (P a ' )l t ) to (P^) 

Before proving our convergence result, we first deduce some simple inequalities 
that will be required. We recall the following well-known results concerning the 
intcrpolant TTh- 

\\(I-nh)4>\\wi.°°iK h )<Ch\tl>\ W 2,<~ (Kh) V0 e [W 2 >°°(K k )] d s x d , 

k = l,...,N K ; (6.1a) 

\\(I-ir h )[ X : 4>]\\mv) < Ch 2 \\Vx\\l^v) ||V0|U~ (2 ,) 

< Cfc Hxlli^z.) ||V^|Uo- ( „) V%, 6 Sft. (6.1b) 

We note for any £ G R>o that 

[^[x:0]] (x) < \ [tth [Cllxll 2 + CW]]0*0 

VxeK k , Vx,0G [C(lQ)} dxd , k=l,...,N K . (6.2) 



i 



i 
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Combining (15TT2")) . (|1.4bl) and (j2l| . we have for all 4> G PZ5 and for all ip G S£ 
that 

' ||^[^/? L (0)]|| 2 < / 7r h [\\tPf3 L m 2 } < f 7r ft [||^|| 2 ||/3 L (0)|| 2 ] 

J D J D 

<dL 2 [ ^[||-0|| 2 ]. (6.3) 

J D 

We require also the L 2 projector IZh '■ V — > V\ defined by 

I (v- K h v)w = VweVi. (6.4) 
Jv 

In addition, we require Vh '■ S — » defined by 

/ 7r h [PfcX:^]= / x:0 V0 G S£ . (6.5) 
Jv Jv 

It is easily deduced for p = 1, . . . , iVp and i, J = 1, . . . , d that 

[PhXh(P P ) = ~t~~ / [PhXlij v P , (6.6) 
Jv Jv 

where r\ v G is such that T) p (P r ) = 5 pr for p, r — 1, . . . , ./Vp. It follows from (|6.5[) 
and (I5.12p with </> = T^x, in both cases, that 

I \\V hX \\ 2 < I M\\V h xf]< I Wxf V X e [L 2 (V)] dx d . (6.7) 
Jv Jv Jv 

We shall assume from now on that T> is convex and that the family {Th}h>o is 
quasi-uniform, i.e. hk > C h, k = 1, . . . , iV/^ It then follows that 

\\K h v\\m(v)<C\\v\\ m(v) WeV, (6.8) 

see Lemma 4.3 in Heywood and RannacherPSEH Similarly, it is easily established 
that 

V X e [H\V)]f d . (6.9) 

Let ([H 1 (V)] d s xd Y be the topological dual of [H 1 (V)} dxd with [L 2 (V)] d s xd being 
the pivot space. Let £ : (\H 1 (V)] d s xd )' -» [H 1 (V)] d s xd be such that £x is the unique 
solution of the Helmholtz problem 

/ [V{£x):-V(t>+{£x)-4>\ = {x,4>)H^v) V0 G [H x {V)]f d , (6.10) 
Jv 

where {-,-)m(v) denotes the duality pairing between ([H 1 (V)] dxd )' and [H 1 (D)]f <d . 
We note that 



X,£x)hhv) = \\£x\\m { v) V X G ([H\V)} d s xd )' , (6.11) 

• _„ / / I f-nMd 

(V) 



and ||£ • is a norm on ([^(D)]^)' 



i 



i 
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Let V be the topological dual of V with the space of weakly divergent free 
functions in \L 2 (V)] d being the pivot space. Let S : V — >• V be such that Sw is the 
unique solution to the Helmholtz-Stokes problem 

/ [V(Sw) : Vv + (Sw) ■ v] = (w,v) v VoeV, (6.12) 

where (■, -)v denotes the duality pairing between V and V. We note that 

(w 7 Sw) v = \\Sw\\ 2 H1{v) VweV, (6.13) 

and \\S ■ || #1(23) is a norm on the reflexive space V. As V is continuously embedded 
in [Hl(V)] d , it follows that [H- 1 (V)] d is continuously embedded in V. 

We recall the following well-known Gagliardo-Nirenberg inequality. Let r € 
[2, oo) if d = 2, and r G [2, 6] if d = 3 and 9 = d(| — j). Then, there exists a 
positive constant C(T>, r, d) such that 



Ivhrp) < C(V,r,mv\\vv,)\\v\\ Hi<D) Vr? e H\V) 



(6.14) 



We recall a lso the follo wing co mpactness result, see e.g. Theorem 2.1 on pl84 in 
Teman JTem84| an d Simon jSim87] Let y Q ^ y anc } y 1 be Banach spaces, y h i = 0,1, 
reflexive, with a compact embedding 3^0^*3^ and a continuous embedding y y%. 
Then, for //j > 1, i = 0, 1, the following embedding is compact : 



{ V €L^(0,T;y ) : ^ G £^(0, T; 3>i) } L»°(Q, T; y) 



(6.15) 



Theorem 6.1. Under the assumptions of Theorem \5.2[ there exists a solution 

{( u a'h ' "'a'DlrSi 6 [Yfc x <$h pdI* °f (P'a h *) suc -h that, in addition to the bounds 
\5.5J$ and \5. 55]) , the following bounds hold: 



=o,...,jv. 



T JT> 



E 



V 



At»a||V<tll 2 



TThl 



a.h 



n = l 



whe 



At n 



#6(2,4) i/d 




and 



i? = 3 if d : 



3.16a) 



h 1 (u) 
(6.16b) 



(6.17) 



Proof. Existence and the bounds (|5.54j) and (|5.55j) were proved in Theorem 15.21 
On choosing <fr = <r a '2 in the version of (|5.50b| dependent on L, it follows from 
(|3~T4|) . (j6T2|l . (|6T3|) . (fS35]) and ([53T]) on applying a Youngs' inequality that 



7T/l 



a.h 



L.n— 1 II 2 l 
r ah II I 



I 



I 
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At n a 

1 

< - 



|Ver 



V 



L,nxx2 
a,h II 



2Wi 



r 1 1 L,n\\2' 

nh[\\<r a \ h \\ ! 



v 



^[ii^r 1 !! 2 : 



2M n f uvu^ii K[<^ l (<'D]II 

JV 



At n 



v 



L , n — 1 1 1 ii »-7 L , n 
\ U a,h \\W V(T a .h\ 



< 



P 



_ r \\L,n-l\\2 i 
^ilka.h II - 



Ai„a 



E E 

t m— 1 p— 1 

■ V—r i: xx' I At 

Ver 



At n d\D\ 
2Wi 



a,?JI 



4Wi 



Ai„C(L) 



IVtt 



L,nii2 
a,h I 



(6.18) 



Hence, summing (|6.18p from n = l,...,m for m = l,...,JVx yields, on noting 
(|5.55[) . the desired result (|6.16a|) . 

n _L,n—l \ 

E V\ in the version of (|5.50aj) depen- 



On choosing w = TZh 



S 



At„ 



dent on L yields, on noting (I6.4[) . (|6.13p . (|6.8|) and Sobolev embedding, that 

2 



Re 



f L.n L,n—1 " 

"»,(■ - u a ,h 



At n 



Re 



v 



L,n L,n — 1 



' L.n L.n — I s 

u a , h ^ u a .h 



At n 



L,n-1 ' 



S 



Re 
~2~ 

Re 
~2~ 



At,, 



u 



r L,n—1 T7\ 
( U a,h ' V ) 



L,n-1 ' 



S 



U a,h - U a,h 



At n 



L,n 



L,n-1 N 



L,n-1 1 



At r 

L,ni|2 



L,ra-li 



It 



-L,ni| ii 2 
a,h II lli 2 (X>) 



IVtt 



L,ntx i|2 



uriiH- 



i(D)J 



(6.19) 



for any # > if d = 2 and for 6 = | if d = 3. Applying the Cauchy-Schwarz and the 
algebraic-geometric mean inequalities, in conjunction with (|6.14[) and the Poincare 
inequality (|1.8[) yields that 



L,n-l| 



, L t n ll l|2 



l W a,'/i"ll lli 2 (r>) - W U a.h 



L,n— 1 ■ ■ 2 



L 4 (2?) H U 'a,/illL 4 (P) 



L,n||2 



<2- E Kniv, 



i 



i 
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n 

<C £ [ll<ri!t^)l|V<rili 2 ( C) ]- (6-20) 



m—n—1 

Similarly, we have for any 6 £ (0, 1), if d = 2, that 



i L.n — In iitt -ZL- ,n 1 1 n2 

Im-A II I V« 1 '- 



"a,h 



^ II L,n — li\2 i|Y7 L,n\\2 

a,h\\ \\L^+f>{V) - \\ U a,h II II VM a,/illL 2 (X>) 



2(1-0) 

<c[[<T ln 1+8 



n 2(l+3g) 

m—n—1 



(6.21a) 



and if d = 3, 

I II L,n-1, 



l x that 

2 

l!(x>) 



IVtt^H II 2 e 



i,n-li|2 
L,n-1 



l o$ llI 2 (X>) 



< cfiKr'iiL-w E n Vw " 



h llL 2 (©)- 



(6.21b) 



m—n — 1 



On taking the J- power of both sides of (|6.19[) . recall (|6 . 1 T[) . multiplying by At„, 
summing from n = 1,...,N T and noting (|6.20|) . (|6.21al) with 9 = |f| i? = 

(Ifi Ol kl\ (IK 07H /1Q 1flnl (IK KKIl /IK Ofll\ /lO /111 ,™1,J„ +U n f 

(i+e) 



a fc*ffl , dOTEb . ([537)1 . (I3~10a . (15351) . (15^91) and (El yields that 



n=l 



L,n 



L,n-1 1 



5 



Atr. 



< CL 2 + C(T) 



N T 



E A< « \\vu L ^\\h {v) + \\r\\ 2 H - 



l(T>) 



1 + max 

n=0,...,A r ' J 



(iKSiii^) 



Nt 



<C(L,T); 

and hence the first bound in (I6.16b[) 
On choosing <fi = Vh 



E A *«n v <tii^w 

,n=0 

(6.22) 



At„ 



G in the version of (|5.50bl) depen- 



dent on L yields, on noting ([63]) . (|6TTD|) . (j6~2l) . (|6~7|) . dHJ), and ([53T]) . that 



L,n 



L.n-l ' 



X> 
1 

Wi 



L,n 



m(r>) 

L.n— I ' 



At K 



' L.n L.n— 1 " 



a.h 



<r \ - er 



L,n-1 ' 



At r , 



I 



I 
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2 / Vu L a ? : n h 



d d 



' L.n L,n—1 " 



At n 



' L.n L,n-1' 
a a.\h ~ a a,h 

At„. 



1 W A £,p( <r a, , fc) 



m=l p=l 



9x E 



/ L.n 



L,n-1 ' 



At r , 



•COL) 



1 + allVcr 



L,n\\1 



L.n II 2 



L,n-ln2 



(6.23) 



'a,ftllL 2 (£>) T II v "a,hllL2(-D) ' ir» a ,h llL 2 (D) 

Multiplying (|6.23[) by Ai„ , summing from n = 1 , . . . , Nt and noting (|5.55l) and 
(|6.16aj) yields the second bound in (I6.16b[) . □ 

6.1. Convergence of the discrete solutions 

First we note the following result. 

Lemma 6.1. For k = 1, . . . , Nk, it follows that 



JK k JK k 



(6.24) 



Proof. We recall the well-known result about equivalence of norms 

i dxd . (6.25) 



1 



7lW<ll*~ sup Uv\\<\\4>\\ V0g 

0,2 vGR d , H| = l 



We recall also that if G 



l dx p d Dl then 



z T 4>z > 



-l it — 1 1 



Vz G 



(6.26) 



that is, H^" 1 !!^ 1 is the smallest eigenvalue of </>. For x S S], PD , on adopting the 
notation in the proof of Lemma 15.11 we have that 

d 

X (x) = J2x(Pj)Vj(x) VxeK k , k = l,...,N K . (6.27) 

3=0 

Then for v G R d , with ||»|| = 1, it follows from ifO?) and flpgj) that 

llx _1 (^) « II > « T xf x (a0 « = CxT 1 ^) «) T x(aO (x -1 ^) v) 

d 

Ellx" 1 (i , *)[[^»7*(*) 



> 



l*- 1 ^)^ 2 Va;etf fe ; (6.28) 



where we have noted that 77^" (a;) > 0, for all x G K k , and x{Pf) G k spl>' J 
0, . . . , d. The bound (gTSS]), on noting (|6T23|l . yields that 

Ilx -1 0»0ll2< [Mllx" 1 !^ 1 ]]^)]" 1 V^Giffc, fc = l,...,iV K , 
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V X G Sl PD . (6.29) 
Hence it follows from (RT2"5|) . $578^ and (JHUI) with n h [ \\x~ 1 \\ ] that 

-i i 



"JX / IIX 

a 2 JJf t 



<f brhiWx-% 1 ]] ^I^IIKIIIX^INIU-C^) 

JK k 

k = l,...,N K , VxeSl,PD, (6-30) 



< C 



Kh[\\X 



□ 



and hence the desired result (|6.24[) . 

We note from (fl~2b|) , $Tfy and that 

II^HII^)]- 1 !! V0GRf^. (6.31) 
Therefore (1535)) . (j6~T6al V). (j5~^9")) . ([532]) . (j6~3T|) . (|63D"|) and (j5~BT)al -c) yield that 



SU P ll M Q,h IIl 2 (-d) 
te(o,T) 



IVti; 



L,At(,±) N 2 



\L 2 (V) 







L, At, +1-1 



i,At,+ 
l a,h 



dt 



L,At,-\\2 



II L,At(,±),|2 
SU P W a ,h IIl 2 (1>) 

*e(o,T) 



IVcr 



L,At(,±)||2 



L.At, 
r i 



du 



L,At 



Of 



lL 2 (X>) 



£■ 



A(t) 



L,At,-|,2 
a,h Hl 2 (X>) 



dt<C , (6.32a) 



L.At 



dt 



A(t) 
d< < C(i,T) 



dt<C(L), (6.32b) 
(6.32c) 



where $ is as defined in (|6.17[) 

We are now in a position to prove the following convergence result for (Pf\ *) 

Theorem 6.2. There exists a subsequence of {(w^'^*, <r^'^*)}h>o.At>Oi ft 71 ^ 
functions u% G L°°(0, T; [Z 2 (2>))] d ) n L 2 (0, T; V) n W 1 ^ (0, T; V) and cr^ g 
L~(0 > T;[i 2 (©)]|^)nL 2 (0 > T;[^ 1 (25)]|?,i) flff^O,! 1 ; ([iJ 1 (2?)]^ xd )') sucft tftat, 
as h, At + , 

weafc* in L°°(0,T; [L 2 (D)] d ), 



L,At(,±) v L 

L,At(,±) v L 
M a ft "> U a 



du 



L,At 
a, ft 



Si 

L,At(,±) 



5 



at 



and 



L.At(,±) L 
r a ,h -> a a 



weakly in L 2 (0,T; [H 1 (V)] d ), 

4 

weakly in L~° (0, T; V), 
strong in L 2 (0,T; [i r (2?)] d ), 

weafc* in L°°(0, T; [L 2 (2?)] dxrf ), (6.34a) 



(6.33a) 
(6.33b) 

(6.33c) 
(6.33d) 



i 



i 
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£,At(,±) L 



L,At 



dt 

L,At(,±) 



<9i 



^[/3 L (<f G±) )]^/? L (^) 
A£, P (<'f (,±) )^/3V£)<W 



wea% in L 2 (0, T; [i7 1 (2?)] dxd ), (6.34b) 

weafc/y m L 2 (0,T; [i? 1 (2?)] dxd ), (6.34c) 

strong in L 2 (0, T; [L r (V)] dxd ), (6.34d) 

strongly in L 2 (0, T; [L 2 (£>)] rfxd ), (6.34e) 
strong in L 2 (0,T; [L 2 (V)] dxd ), 

m,p = l,...,d, (6.34f) 



where d is defined by l[6.17\ ) and r G [1, oo) if d = 2 and r £ [1,6) if d = 3. 
Furthermore, (u^,(T^) solve the following problem: 

(P^) Find u L a G L°°(0,T; [£ 2 (£>)]< i ) PI L 2 (0,T; V) n W^ 1 '^ (0, T; V) and cT L a G 
L°°(0,r[L 2 (D)]f P d D )nL 2 (0,r; [77 1 (^)]5PD)ni/ 1 (0,T;([i/ 1 (^)]5 >< ' i )') «*cA i/iai 



Re 



dt 



dt 



f f [(1 - e) Vu£ : + Re ■ V)u L a ] ■ v] , 
= [{f,v) Hl{v) dt-^J T Jj^):Vvdt 



JX> 

W G LT^(0,T;V), 



(6.35a) 



m(v) 



f f [(^■V)^(^)]:^ a V^:: V0] 



'0 JD 



2(Vu£)0 £ («t£)- — {a L a -I) 



4>dt 



G i 2 (0,T; [ff 1 (X>)]| xd ) 



(6.35b) 



and 



lim / («*(i,x)-tt°(x)) • t7 = 

Mv G H := {tt> G [L 2 (2?)] d : divw = in V} , 
lim / (<r L a (t, x) a°(x)) : X = V* 6 [L 2 (V)] d % . (6.35c) 



Proof. The results ()6.33al -c) follow immediately from the bounds (|6.32a[ c) on not- 
ing the notation (|5.60al -c). The denseness of {J h>Q Q,\ in L 2 (D) and (|5.1dj) yield 
that u L a G L 2 (0,T;V). The strong convergence result (|6.33d|) for u A ' follows 
immediately from (|6.33al -c) and ([6351) with ,uq = 2, = 4/i9, >b = [-ff 1 ( :z? )] d , 
= V with norm ||<S • ||ijim) and ^ = [L r iV)] d for the stated values of $ and r. 
Here we note that y is a Banach space and X> * = 0, 1, are reflexive Banach spaces 
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with [L r (T>)] d continuously embedded in V', as [H~ 1 (T>)] d is continuously embed- 
ded in V', and [H 1 (T>)] d compactly embedded in [L r (T>)] d for the stated values of 
r. We now prove (l6~33d)) for u^f First we obtain from the bound on the last 
term on the left-hand side of (I6.32a|) and (|5.61[) that 

ll u a,'h ' ~ u a,h t ' ± lli=(o,r,L=(a?)) ^ c At - ( 6 - 36 ) 
Second, we note from Sobolev embedding that, for all r\ G L 2 (0, T; 7J 1 (2?)), 

h\\Li{Q,T;Lr{V)) < \\v\\l^(0,T;L^(V)) IMIz,=(0,T;£»(X>)) 

< c \\v\\l"(o,T;L"(t>)) hWl^o, T-HHV)) ( 6 - 37 ) 

for all r G [2,s), with any s G (2,oo) if d = 2 or any s G (2,6] if rf = 3, and 
= [2 (s - r)]/[r (s - 2)] G (0, 1]. Hence, combining $6M§ . (f6737|) and (|633d|) for 

"a'h* yields (|6.33dj) for u^*' . 

Similarly, the results (|6.34al -c) follow immediately from (|6.32b| c). The strong 
convergence result (|6.34d|) for it^j' follows immediately from (|6.34al -c). (|6. 1 1|) and 
(lQ5l) with fi a = (Mi =2,y a = [H 1 (V)] dxd , = [H- 1 (V)] dxd and y = [L r (V)] dxd 
for the stated values r. Similarly to ()6.36|) . the second bound in (|6.32b|) then yields 
that (|Q4d)l holds for cr^f t(,±) . 

Since v^'h G L 2 (0,T; S\ PD ), it follows that er^ is symmetric non-negative 
definite a.e. in T>t- We now establish that <x^ is symmetric positive definite a.e. 
in T>t- Assume that <x^ is not symmetric positive definite a.e. in V\ c T>t- Let 
v e L 2 (0,T; [L 2 (V)] d ) be such that a^v = with ||«|| = 1 a.e. in V% and v = 
a.e. in £> \ "Dj,. We then have from (|6.32aj) that 

<C^£ J^ L a f>+,(vv T )dtY . (6.38) 

Hence it follows from (|638]) and (|6.34d|) that |2>£| = 0. 

Finally, the desired results (|6.34cl f) follow immediately from (j5.64[) the second 
bound in (jQ2bll . (l2~T71) . (jQ4dll and the fact that cr^ G L°°(0,T; [L 2 (£>)]^). 

It remains to prove that (u^,<x^) solves (P^). It follows from f|5 .3(1 . (|6.32al -c). 
(|6.33aJ -d). (|6.34e|) . (|3.10bj> . (f6TT2|) and ([3~T3l> that we may pass to the limit, h, At -> 
0+, in the L-dependent version of (|5.62j) to obtain that (u„,(T„) satisfy (|6.35a[) . It 
also follows from (|5.28a|l . ([53|1 . (|6.33cl d) and as V is dense in H that tt ^(0, ■) = u°(-) 
in the required sense; see (|6.35cl) and Lemma 1.4 on pl79 in Temam! em " I 

It follows from (|Q3aT -f). (|6.33bl d), ([BTTD]) . (|Q2lil -c). ([BTTal b). IfLla]) and as 
ix^ G L 2 (0, T; V) that we may pass to the limit /i, At — > 0+ in the L-depcndent ver- 
sion of (15^31 with x =n h (f>to obtain (l6~35bl) for any d> G C£°(0, T; [C°°(25)]g x<< ). 
For example, in order to pass to the limit on the first term in the L-dependent 



i 



i 
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version of (|5.63|) , we note that 

f-T 



JT> 



d(T ah 1 



L,At. 



Of 



Wi 



T 

JV 



dt Wi <Tq <' 1 



1 



TT h (j) dt 

TT h (p + (I - TT h ) 



L,At 



d<t> 



dt 



Wi 



JT> 



L,At,+ 



dt. 



dt 



.39) 



The desired result (|6.35b|) then follows from noting that C$°(Q,T; [C°°{V)]g) is 
dense in L 2 (0,T; [H 1 {V)]f- d ). Finally, it follows from (IE2lbl) . (|53Icl dL (IBTal b) 
and (|5. 12[) that er^ (0, ■) = <r°(-) in the required sense; see (|6.35c|) and Lemma 1.4 

□ 



pl79 in Temaml^SmHl 



Remark 6.1. It follows from (|6.32a| b). (|6.33al b) and (|6.34al b) that 



SUP 

te(o,T) 



|VuSlli=.(p)*<C i 



(6.40a) 



sup II^Hl^+a/ |]V^||| 3(P) df<C(i). (6.40b) 
te(o,T) Jo 

Hence, although we have introduced a cut-off L > 1 to certain terms, and added 
diffusion with a positive coefficient a in the stress equation compared to the stan- 
dard Oldroyd-B model; the bound (|6.40a|) on the velocity is independent of the 
parameters L and a, where (u^,cr^) solves (P^), (|6.35al -c). 

7. Convergence of (P^' h ) to (P a ) in the case d — 2 

First, we recall the discrete Gronwall inequality: 
(r°) 2 + (s o r<( q °) 2 , 



m— 1 77i 



m > 1 



TO > 1 . 



(7.1) 



n=0 



n=0 



Theorem 7.1. Under the assumptions of Theorem \5.2\ there exists a solution 
{(«S,h.<h)}»=i e [VJ x S^]^ o/ fifj ^ *e bounds f£5fl and 

nrb%) hold. 

If d = 2, a < and At < C+(( : ~ 1 ) a 1+ ^ h 2 , for a £ > 0, then the following 
bounds hold: 



max / ^MK.J 2 : 

n=0,...,N T J v 



N T 

E 

n=l 



Ai„a||V<.,J| 2 



•^[iKfc-^ii 3 : 



i 



i 
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N T 



E A ^ 

n=l 



s 



At r 



<C{a-\T)- (7.2) 

ffl(Z>) 



where $ € (2,4). 



Proof. Existence and the bounds (|5.54j) and (15.551) were proved in Theorem 15.21 

On choosing <fi = cr 1 ^ h in the //-independent version of (|5.50bj) . it follows from 
(|3.14p and on applying a Youngs' inequality for any ( > that 

\f **[ IKJI 2 + \\<,h - <;h\? \ + &t n aj^ ii v<,f + ^ jf 7r h [ ikji 2 ] 



< ^ jf «X[ IKi 1 II 2 ] + + 2At„ jf V< h : , h [{K. h f] 

* d d Ocr n 
+ At n £Ka L VpKO : 

^ rn— 1 p— 1 ^ 

< 5/ ^[IKr^fl + c^ [i + ||v<,,iu 2( p) \\Kh[{<, h f]\\L*(p)] 



.n-li 



+ C At„ ||tt™ ^ || 2(2+Q 

ll A m,p(° r a,/ l )llL2+C(I5) || V<T™ fe || L 2(p) . (7.3) 
L C (U) 

It follows from (|5TT2j) . (ll~4b|) . dHJ) and (|oTT4]) . as d = 2, that 

\^[« h ) 2 }\\h { v) = J v Uh\{< lh f] f < J v m \\(<, h n 2 }<f v Kh[\\cr n a , h \\ 4 } 

N K . N K 

= E / ^[IKhll 4 ]<El^llKhllt-(x fc ) 



< C J2 \Kk\ (\Kk\~ l IKJIiH^)) 4 < C E H<J V*) 
fc=l fc=l 

= C ll CT 2,/Jlz,4(z>) < C Ika^lli^cD) ■ (7-4) 

Similarly, it follows from the 5- independent versions of (|5.22l) . (IS^lal b - ). recall Re- 
mark O (|5T2l3j) . (g3) and (joTTl) that for all C > 

N K N k 

\\^A< h )\\l + 4 HV} < E 1**1 IIA m , P «,,)|iit c (Kfe) < c E 1^*1 IKfcll^*) 
fe=i fe=i 

= c iks,hll^c (D ) < c (0 IKhlli» (JI ) IK Ji^ ■ (7-5) 

In addition, we note from (foTT4]) , (fT5]l and (|535j) that for all C > 



L <• (X>J 



(7.6) 



I 



I 
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Combining (fOjl . (fTi]) . ([731) and ((776]), and on noting (|5~T2|) and that a < 
yields on applying a Young's inequality that for all ( > 

At, 



^[IKJI 2 



,n-l 1 1 2 - 



Ai„a 



IVo-" J| 2 



< / 7r h [||<^|| 2 ]+CAt„ 
+ C(r 1 )At„a-( 1+ « 



2Wi 



l^ M S,/illz 2 (T>) 



"7«[ ll CT a,ftll J 



(7.7) 

Hence, summing (I7.7[) from n = 1, . . . , to for to = 1, . . . , iVy yields, for any £ > 
that 



™- r ii ~ m ii 2 " 

^1 ll^/JI J 



m ~ 

E / 



«||V<J 



1 



2Wi 



E 



n=l 



„. r II _n n — 1 1|2 " 

^[ll^a.h II . 



< / Kh[\K h \n+c 



E A< « 



I ^NKJ 2 ]. (7.8) 



E n v <^i 

_k— n — 1 

Applying the discrete Gronwall inequality (|7.1j) to fjT. 8[) . and noting (|5.27[) . (|5.13|) . 
(15^91) . (l5~55jl . (1531) and that At < C^C 1 ) a 1+<: h 2 , for a £ > where C^C -1 ) is 
sufficiently small, yields the first three bounds in (|7.2j) . 

Similarly to (|6.19[) . on choosing w = TZh 

independent version of (|5.50ap yields, on noting 
embedding, that 

2 



S 



At„ 



G Y\ in the L- 



(j633l) . (EHI) and Sobolev 



Re 



U a,h U 'a,h 



Re 



Ai, 



s 



"a. ft, "a, ft 



At r , 



(7.9) 



< c[\wi h \\i 2(v) + iivu^iil,^) + 1| k^h ||<J 
+ 11 IK^IIIIvk,,!! \\l +a(v) + \\n 2 H - HV) ] 

for any > as d = 2. On taking the 1 power of both sides of (|7.9I) . multiplying by 
At„, summing from n = 1, . . . , Nt and noting the L-independent versions of (|6.20[) 
and (fOTal) with 9 = (?? - 2)/(6 - tf), (|Q7) . (ETTOal) . (15351) . (l5~29l) . ([5~T2)l and the 
first bound in (|7.2p yields the last bound in (|7.2j) . □ 



It follows from (|535|) . ([7^]) . (j5~^9"|) . ([5~T2]) . (jOO)) and (j53Dal -c) that 
T ||V« At „ ( ' ±) ||i 2(I , ) dt 



II At(,±) 1,2 
SU P IK.ft WtfCD) 

te(o,T) 



June 23, 2010 6:7 jwbsblarxiv2 



Existence and Approximation of a (Regularized) Oldroyd-B model 49 



u 



At,+ 



U 



At,— 1.2 

aji \\L 2 (T>) 



A(t) 



dt<C (7.10a) 



and 



I At(,±),|2 
SU P W^a.h h'fp) 

te(o,T) 



allV^l^" 2 



r A£,+ 

a,h 



L 2 (V) 



r At,~,,2 
T aji \\L 2 (T>) 



A(t) 



+ 



5 



r)n At 
ou aJi 



dt 



dt < C(a-\T), 



(7.10b) 



where & G (2, 4). 

We note that we have no control on the time derivative of (T^\ in (|7.10bp . 



This is because if we choose <p = Vh 



G in the L-independent 



version of (|5.50b[) . the terms involving u™ h , m = n — 1 and m = n, cannot now be 
controlled in the absence of the cut-off on <x™ h . We are now in a position to prove 
the following convergence result for (P A t h ). The key difference between the following 
theorem and Theorem 16.21 for (Pfj'i^*) is that no control on the time derivative of 
<r^ in (17.10b)) implies no strong convergence for 0"^ /» • 

Theorem 7.2. Let all the assumptions of Theorem \7.1\ hold. Then there exists a 



subsequence of o - „ i /,)}h>o,At>o, and functions u a e L°°(0, T; [L 2 (£>))] 2 ) n 

L a (0,T;V)nW 1 '*(0iT;V ; ) and cr & G L°°(0,T; [L 2 (2?)] 2 s ^)nL 2 (0, T; [H 1 ^)]^) 
suc/i i/iai 7 as /i, At — > + , 



At(,±) 



At(,±) . 



du 



dt 

At(,± 
l a,h 



h c du. 

► S 



dt 



and 



At(,±) 



->■ cr n 



At(,±) 



A 



m,p 



( Ai(,±)s 



iii'P 



weak* in L°°(0, T; [£ 2 (P)] 2 ), (7.11a) 

weafcZy in i 2 (0,T; [^(D)] 2 ), (7.11b) 

wea% m L*(0,T; V), (7.11c) 

strongly in L 2 (0,T; [L r (2?)] 2 ), (7.11d) 

weafc* m L°°(0,T; [L 2 (D)] 2x2 ), (7.12a) 

wea% in L 2 (0, T; [iJ 1 (2?)] 2x2 ) I (7.12b) 
weafcfy in i 2 (0,T; [L 2 (£>)] 2><2 ), 

m,p=l, 2, (7.12c) 



where G (2, 4) and r G [1, oo). 

Furthermore, (u a ,cr a ) solve the following problem: 



I 



I 
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(P a ) Find u a G L°°(0,T; [L 2 {V)} 2 ) nL 2 (0,T;V) PI W 1, $(Q,T; V') and er a G 
L°°(0,T; [L 2 (V)] 2x2 D ) n L 2 (0,T; [H l {V)f s %) such that 



Re 



"at 



dt 



[(1 - e) Vu a : Vv + Re [(tt a • V)tt Q ] • «] dt 



J o (f,v) H i {v) dt 



e 

Wi 



Vvdt 



v 



Vw 6 L— (0,T;V), 
(7.13a) 



#0 
at 



dt - 



4> 



[ [ [{u a -V)(r a : 4> + aV(r a :: Vd>] dt 

Jo Jv 



2{Vu a )(T a - —{<T a I 

Wi 



lim 



J + Jv 



V0 G L 2 {Q,T;[H l (V)fg 
(u a (t, x) — u [x)) ■ v = 



4>dt 

-T,T;[L 2 (V)] 2 s x2 ); (7.13b) 
V«eH. (7.13c) 



Proof. The results (|7.11af c0 and (|7.12al b) follow immediately from the bounds 
^.lOal b) , as in the proof of Theorcm l6.2l Similarly, the proof of positive definiteness 
of cr a follows as in Theorem 16.21 that is, (|6.38l) and the weak convergence (|7.12a|) 
is adequate for this. The result (|7.12c|t follows from (|7.12ajl . ([5lM| and (|7.10b[) and 
the fact that er Q G L°°(0,T; [L 2 {V)] 2 s x2 D ) . 

It follows from ([O]), (TTIUal b). (fHTal -dV (|7l2a|l . (l3~10bl) . (IcTT^ and (13TT31) 
that we may pass to the limit, h, At — > 0+, in the L-independent version of (|5.62|) 
to obtain that (u a ,cr a ) satisfy (|7. 13a|) . It also follows from (I5.28a|) . (|5.3p . (JTTTcld) 



and as V is dense in H that tt ^fO, •) = u°(-) in the required sense; see (I7.13cp and 
Lemma 1.4 on pl79 in Temamff em ^ 4 l 

It follows from (j7T2U -c). (|7.11dD . (jcTTUl) . (ITTTUal b). (IrTTal b). (H~4aj) and as u a G 
L 2 (0, T; V) that we may pass to the limit h, At —> 0+ in the L- independent version 
of (15^31) with x = n h 4> to obtain (l7TT3bl) for any G C^°(-T,T; [C 00 ^)]^ 2 ). 
For example, in order to pass to the limit on the first and third terms in the L- 
independent version of (|5 .63[) , we note that 
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dt 



dt 
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rr At ■ or 
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dt. 



0] (0, •) 
(7.14a) 
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/ vr At.+ r At,+ M 

/ VU ah :TT h [<T h 7T h( j)} 

Jv 

= I v<,; + :(^-/)[<r^</>] 

Jv 

- I {((V** <, + • ((-^)div<r)} ; (7.14b) 

where ((Vir h <j))u^\ + )(t,x) G K 2x2 with 

[(V7r/,^)u££% = ELi [<h + ]fc- The desired result then fol- 

lows from noting that Cg°(-T,T; [C°°(V)] 2x2 ) is dense in W^'\0,T; {H^V)} 2 * 2 ^ 

We have the analogue of Remark 16.11 
Remark 7.1. It follows from (|7.10a[ b). (|7.11al b) and (|7.12a( b) that 

sup Hitall^m) + / ||Vu Q ||i2 (2J) dt<C, (7.15a) 
te(o,T) Jo 

sup [|o- a ||| 2(P) +a f WVcTcWl^dtKCia-^T). (7.15b) 
te(o,r) Jo 

Hence, although we have introduced diffusion with a positive coefficient a into 
the stress equation (|7.13b[) compared to the standard Oldroyd-B model; the bound 
(|7.15al) on the velocity u a is independent of the parameter a, where (u a ,cr a ) solves 

(p ), dzna-c). 
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